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Abstract 

^ • Originally motivated by a stability problem in Fluid Mechanics, we study the spectral 

I and pseudospectral properties of the differential operator iJ^ = —d^ + x'^ + ie~'^f{x) on 

l/^ ■ L'^(M), where / is a real- valued function and e > a small parameter. We define S(e) as the 

I infimum of the real part of the spectrum of H^, and 5'(e)~^ as the supremum of the norm of 

^\ • the resolvent of -ffc along the imaginary axis. Under appropriate conditions on /, we show 

' that both quantities S(e), ^'(e) go to infinity as e ^ 0, and we give precise estimates of the 

OO i growth rate of ^'(e). We also provide an example where S(e) ^ ^(e) if e is small. Our main 

results are established using variational "hypocoercive" methods, localization techniques and 
semiclassical subelliptic estimates. 



c5 : 1 Introduction 

In many evolution equations arising in Mathematical Physics, one encounters the situation where 
the generator of the evolution can be written as a sum of a dissipative and a conservative operator 
which do not commute with each other. In such a case the conservative term can affect and 
sometimes enhance the dissipative effects or the regularizing properties of the whole system. For 
instance, if the system has a globally attracting equilibrium, the rate of convergence towards 
this steady state can strongly depend on the nature and the size of the conservative terms. 
Typical examples illustrating such an interplay between diffusion and transport are the kinetic 
Fokker-Planck equation [13], and the Boltzmann equation [5]; see also [21] for a comprehensive 
study of these phenomena at a more abstract level. 

In this paper we study a simple linear system which fits into this general framework. Given 
a small parameter e > and a smooth, bounded function / : M — > M, we consider the differential 
operator 

= -d'', + x^+'-f{x) , x£R, (1.1) 
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acting on the Hilbert space X = ^^(M), with domain V = {u G H'^{M);x^u E ^^(M)}. 
Clearly iJ^ is a bounded, skew-symmetric perturbation of the harmonic oscillator Hao = —d'^+x^. 
Our goal is to compute the decay rate in time of the solutions to the evolution equation 

All 

— = -H,u , u{0) = no e X . (1.2) 

As we shall see, the solutions to ()1.2p decay to zero at least like e~*' as t — > +00, but the actual 
convergence rate strongly depends on the value of e and the detailed properties of /, due to the 
interaction between the symmetric (dissipative) and the skew-symmetric (conservative) part of 
the generator —H^. 

Our initial motivation for this study is a specific problem in Fluid Mechanics which we now 
briefly describe. As is explained in [6l[7], to investigate the long-time behavior of solutions to the 
two-dimensional Navier-Stokes equation, it is convenient to use the the vorticity formulation. In 
self-similar variables, the system reads: 

+ u ■ Vuj = Auj + -X ■ Vuj + uj , 2; G , i > , (1.3) 

where uj{x,t) G M is the vorticity distribution and u{x,t) G M? is the divergence- free velocity field 
obtained from to via the Biot-Savart law. Equation (jl.3p has a family of stationary solutions, 
called Oseen vortices, of the form uj = aG where G{x) = (47r)~-'^e~l^l and a G M is a free 
parameter (the circulation Reynolds number). It turns out that the linearization of (jl.Sp at 
aG has the same form as (|1.2p . in the sense that the generator can be written as a difference 
C — aA, where £ is a self-adjoint operator in the weighted space L^(M^, da;) and A is a 
skew-symmetric perturbation. The analogy goes even further if we conjugate C — aA with the 
Gaussian weig ht and if we neglect a nonlocal, lower-order term in the perturbation A. The 
linearized operator then becomes 

^„ = -A + ^ - i + af{x)de , X G R^ (1.4) 
lb 2 

where = xid2 — X2di and /(x) = (27r|xp)~^(l — e"'^''^/^). The operator in (jl.ip is a one- 
dimensional analog of Ha, and the limit e ^ corresponds to the fast rotation limit a — > 00. 
Remark that, in this particular example, the function / has a unique critical point located at 
the origin, and decreases to zero like as |x| 00. 

The aim of this paper is to study the spectral and pseudospectral properties of the linear 
operator in the limit e ^ 0. Besides the specific motivations explained above, this question 
has its own interest from a mathematical point of view, and turns out to be relatively complex. 
We have to deal with a non-self-adjoint problem of (almost) semiclassical type which exhibits a 
competition between various microlocal models at different scales, depending on the structure 
of the critical points and the decay rate at infinity of the function /. In particular, unlike in the 
self- adjoint case, or even in some non-self- adjoint problems such as the kinetic Fokker-Planck 
operator (see [TOl EHl E] ) , the pseudospectral estimates are not monotone with respect to the 
imaginary part of the spectral parameter. Nevertheless, the model (jl.ip is simple enough so 
that the analysis can be pushed quite far, and we believe that our results give a good idea of 
the phenomena that can be expected to occur in more general situations. We also mention that 
the spectral theory of non-self-adjoint operators, especially in the semiclassical limit, is a topic 
of current interest [23l [3l lil [22] . 

We start with a few basic observations concerning the operator H^. As is well known the lim- 
iting operator H^c = —d^ + x'^ is self-adjoint in L^(M) with compact resolvent, and its spectrum 
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is a sequence of simple eigenvalues {A^}„gN, where = 2n + 1. By classical perturbation the- 
ory [T7j, it follows that has a compact resolvent for any e > 0, and that its spectrum is again 
a sequence of (simple) eigenvalues {An(e)}neN5 with Re(A„(e)) — > +00 as n ^ cxd. Moreover, it 
is clear from (jl.ip that the numerical range 

Q{H,) = ^{H,u,u)l2 £C;u£V, ||n||i2 = l} (1.5) 

is contained in the region TZ^ C C defined by 



A G C; Re(A) > 1, elm(A) G /(M)} . (1.6) 



In particular we have A„(e) G TZ^ for all n G N and all e > 0. The operator is therefore 
sectorial, and — 1 is maximal accretive, see Section V-3.10 in |17j . By the Lumer-Phillips 
theorem, it follows that 

||e"*^^|| < e"* , for all t > . 

Here and in the sequel we denote by ||^|| the norm of a bounded linear operator A in L^(M), 
and by Rsp{A) its spectral radius. 

Since we want to compute the decay rate of the solutions to (II. 2|) . we are interested in 
locating the real part of the spectrum of H^. We thus define: 

S(e) = miRe{a{H,)) = minRe(A„(e)) . (1.7) 

neN 

We also introduce a related pseudospectral quantity, which will play an important role in this 
work: ^ 

^(e) = (sup\\{H,-i\)-'\\) . (1.8) 
It is easy to see that S(e) > ^'(e) > 1. Indeed, for any A G M, we have 

^istaxMH.)) = ^.((^.-'^)-'> i UH.-.xn\ < ,i,,(,,;e(g.)) ■ (i-»> 

see e.g. Problem III-6.16 and Theorem V-3.2 in ]17| . Taking the supremum over A G M and 
using the inclusion Q{Hf^) C TZ^, we thus find 

1 1 1 



S(e) - ^'(e) - ASM dist(a,7^e 

which is the desired result. Elementary relations between S(e), 'J'(e) and the norm of the 
semigroup e~*^' can also be deduced from the following general lemma: 

Lemma 1.1 Let A be a maximal accretive operator in a Hilbert space X, with numerical range 
contained in the sector {z G C ; |argz| < ^ — 2a} for some a G (0,^]. Assume that A is 
invertible and let 

S = miRe{a{A)) > , and ^' = fsup ||(A - iA)"^||) 

VAeM / 

Then the following holds: 

i) If there exist C > 1 and /i > such that ||e~*^|| < C e"'^* for all t > 0, then 

^ - ^ ' ^ - l + log(C7) ■ 
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ii) For any /i G (0, S), we have \\e ^"^W < C{A, jj) e ^* for all t > 0, where 

C{A,ij) = — - — (i2N{A,fi) + 2TT] , and N{A,fi) = sup ||(^ - ^ - iA)"^|| . 
vr tan a V / AeK 

in) If moreover fi G (0, the quantity N{A,fi) is not larger than — fj.)^^. 

For completeness, we give a proof of this lemma in Appendix |Al We also mention that a 
similar result holds for some maximal accretive operators with a = 0, in which case the constant 
C{A, n) has a different expression, see [13] or [10] for a detailed study of Fokker-Planck operators 
which enter in this category. In the case when A = H^, we have of course S = S(e), ^ = ^'(e), 
and the angle a satisfies tan(2a) = e||/||^i, see (|1.5p . (|1.6p . Thus a = 0{e) as e — > 0. 

Lemma 1 1.1 1 can be used in particular to illustrate the pseudospectral nature of the quantity 
^'(e) defined in (jl.Sp . For semiclassical operators of the form = p{x, edx), the pseudospectrum 
is usually defined as the set of all z G C such that e^||(Pe — z)""*^ || — > oo as e — > for all G N, 
see e.g. [19]. Instead of defining the pseudospectrum as a subset of C, we find it convenient to 
introduce here a more flexible notion: 

Definition 1.2 Let (we)gg(o,i] ^ family of complex domains, i.e. C C for all e G (0, 1]. 
We say that We meets the pseudospectrum of as e ^ if 

lime^ sup \\{H, - z)-^\\ = +oo , for a// A G N . 

On the contrary, we say that uj^ avoids the pseudospectrum of as e ^ if there exists A G N 
such that 

sup\\{H,- z)-^\\ = 0{e-^) , ase^O. 

If we apply this deflnition in the case when = {z G C ; Re(2;) < /x^}, we arrive at very 
different conclusions depending on whether /x^ <C ^(e) or fj,^ ^ ^(e)- Indeed, using Lemma ll.ll 
and Deflnition ll.21 we obtain the following result, whose proof is again postponed to Appendix lAl 

Lemma 1.3 

i) For any k G (0, 1), the domain {Re(2;) < K^'(e)} avoids the pseudospectrum of as e — > 0. 

ii) If fie ^ + log ^(e) + log(e~"'^)) in the sense that the ratio goes to +oo as e — > 0, then 
the domain {Re(z) < /x^} meets the pseudospectrum of as e ^ 0. 

The main purpose of the present work is to investigate the behavior of both quantities S(e) 
and ^'(e) as e ^ 0. In particular, we shall give precise estimates on ^'(e) under rather general 
assumptions on the function /. Computing S(e) is a more delicate task, and we shall restrict 
ourselves to a specific example which shows that S(e) can be much larger than ^'(e) when e 
is small. This can be expected because H^. is highly non-self-adjoint in this regime (we recall 
that S = ^' for self-adjoint operators). 

If the function / : M ^ M is not identically constant, it is rather straightforward to verify 
that S(e) > ^(e) > 1 for all e > 0, see Lemma |2. II below. In view of Lemma ll.ll we conclude 
that the decay properties of the semigroup e~^^^ are always enhanced by the skew-symmetric 
perturbation ie~^f{x), if / is not constant. It is therefore very natural to ask under which 
conditions on / the decay rate can be made arbitrarily large by choosing e sufficiently small. If 
we use the pseudospectral quantity ^'(e) to measure the decay, the complete answer is given by 
the following proposition, which is proved in Section [2l 



4 



Proposition 1.4 Assume that f € L°°(M) n C'^(M). // all level sets of f have empty interior, 
then ^'(e) ^ oo as e —> 0. In the converse case, ^(e) is uniformly bounded for all e > 0. 

We suspect that Proposition II .41 remains true if we replace ^'(e) by S(e), but at the present time 
there is no proof that S(e) remains bounded if / is constant on some nonempty open interval. 

We now turn to more quantitative versions of Proposition [L4l which specify the growth rate 
of ^(e) or S(e) as e ^ 0. Our first result in this direction will be obtained using a variational 
method recently developed by C. Villani under the name of "hypocoercivity" ^24j. This approach 
applies (in particular) to linear operators of the form L = A* A + i? in a Hilbert space X, where 
B is skew- symmetric. Under certain conditions, it allows to compare the spectral properties of 
L with those of the self-adjoint operator L = A* A + C*C, where C = [A, B]. It is clear that our 
problem fits into this general framework: if we set X = L^(M, C), ^ = dx+x, and B = {i/e)f{x), 
we see that H^ — l = A*A+B. Since C = [A, B] = (i/e)/'(x), the associated self-adjoint operator 
He defined hy — 1 = A* A + C*C has the exphcit form 



= -dl + x^ + \f'{xf , X G M . (1.10) 



The point is that the spectral properties of in the limit e — > are rather easy to study using 
semiclassical techniques, see e.g. Lemma [1 .71 below. Once the spectrum of is known, Villani's 
method allows to deduce useful information on the original operator H^. A great interest of this 
approach is that it can be easily adapted to nonlinear problems [2U [25] . 

Following closely the proof of the "basic hypocoercivity result" in [231 Section 4] , we obtain 
our first main result: 

Theorem 1.5 Assume that f G C^(M) satisfies f",f"' G L°°(M), and that there exist Mi > 
and V G (0, 1/2] such that 

{HeU,u)L2 = [ (|5,u|2+xV|2 + ^/'(x)>|2)dx > ^ , (1.11) 

for all u £ T> and all e G (0, 1]. Then there exists M2 > such that, for all e G (0, 1], 

S e > — , and ^(e > — — • 1.12) 

Remark that we do not suppose in Theorem 1 1 . 5 1 that / is bounded, but only that /" and 
are. Under these more general assumptions, the perturbation B.^ — B^q is relatively bounded 
with respect to -ffoo, but not necessarily relatively compact. To test the conclusion on a concrete 
example, consider the simple case when f{x) = 'j? . Then we have B^ = —d^ + (1 -|- i/e)x^, so 
that S(e) = Re((l + Ve)^/^) = 0{e-^/^) as e ^ 0. On the other hand, B, = -dl + {I + A/e^)x^ , 
hence inf(cr(i?,)) = (1 -h i/e^)'^/'^ so that (fTTTI) is satisfied with 1/ = 1/2. Thus the lower 
bound (jl.l2p for S(e) is optimal, and one can also check that the estimate for ^'(e) is optimal 
except for the logarithmic correction. However, as we shall see, the results given by Theorem 1 1.5 1 
are not always optimal, especially when f'{x) decays algebraically as — > 00. 

Motivated by our original example ()1.4p . we next consider a specific class of nonlinearities 
with a prescribed behavior at infinity. 

Hypothesis 1.6 We assume that f G C^(R,M) has the following properties: 
i) All critical points of f are non- degenerate; i.e., f'{x) = implies f"{x) / 0. 
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ii) There exist positive constants C and k such that, for all x GM. with \x\ > 1, 



X 



< , for £ = 0,1,2,3 . (1.13) 



X 



Loosely speaking, we consider Morse functions which are bounded together with their deriva- 
tives up to third order, and which behave hke \x\~^ as oo. This is not the most general 
class that we can treat with our methods, but it already contains many interesting examples, 
including the original problem (|1.4p . The assumption that f'{x) decays algebraically as |x| oo 
is very important: as we shall see, the decay rate of the semigroup e~*^' is determined not only 
by the structure of the critical points of /, but also by the behavior of / near infinity, and under 
hypothesis (jl.lSp there is a competition between both regimes which is one of the main interests 
of our problem. 

Under Hypothesis 11.61 it is rather straightforward to estimate the lowest eigenvalue of the 
self-adjoint operator iJ^. The result is: 

Lemma 1.7 /// satisfies Hypothesis \ 1.61 there exists M3 > 1 such that, for all e € (0, 1], 



,^ „ < inf crfi^J < ^ , where v = • (1-14) 

For the reader's convenience, we sketch the proof of this result in Appendix [Al Since (jl.l4p 
implies (jl.lip with the same value of u. Theorem 11.51 shows that S(e),^'(e) obey the lower 
bounds (|1.12p with v = {k + 2)~^ if / satisfies Hvpothesis 11.61 It turns out that these estimates 
are not optimal in the present case, neither for $](e) nor for ^'(e). As is discussed in Section [3] 
below, this is not an intrinsic limitation of the variational method, but the proof of a general 
result such as Theorem 11.51 fails to take into account all specificities of /. Our second main 
result gives optimal estimates on ^'(e) under Hvpothesis ll.6t 



Theorem 1.8 /// satisfies Hypothesis \1.6l there exists M4, > 1 such that, for all e G (0, 1], 

< ^(e) < , where v = • (1.15) 



Remark that z7 > z/, so that the lower bound in (jl.lSp is strictly better than in (jl.l2p . In this 
paper, we shall give two different proofs of this lower bound. The first one consists in adapting 
the variational method that we already used in the proof of Theorem 11.51 This approach gives 
the lower bound in ()1.15p with a logarithmic correction, as in (jl.l2p . The second proof, which 
removes the logarithmic correction and also provides the upper bound in (jl.lSp . uses rather 
classical techniques in the analysis of partial differential operators. The idea is first to reduce 
the problem to a bounded domain using a dyadic partition of unity in M, and then to apply 
semiclassical subelliptic techniques as developed, for instance, in [1] and in Chap. 27 of [16] . 

Finally, we discuss the behavior of the spectral quantity S(e). Under Hvpothesis 11.61 Theo- 
rem [L8l already provides a good lower bound on S(e), via the inequality $](e) > ^'(e). In fact, 
if one believes that the analysis of small random perturbations of non-self-adjoint semiclassical 
operators presented in [8| can be carried over to the operator H^, one also expects that the lower 
bound S(e) > ^'(e) is optimal for a generic / in a wide class of functions. However S(e) can be 
much larger than ^'(e) in some particular situations. This is the case for instance when f{x) = x, 
because one can show that S(e) = ©(e'^) and ^(e) = ©(e'^/^), see Section [322] below. To 
conclude this paper, we consider a more interesting example which satisfies Hypothesis 11.61 and 
is inspired by the original problem (jl.4p . 
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Proposition 1.9 Fix k > and assume that 

Then there exists a constant M5 > such that the lowest real part of the spectrum satisfies, for 
alUe (0,1], 

S(e) > — r , where v' = min-f— , \ . (1-17) 

Since v < u' , it follows from (|1.15p and (|1.17p that S(e)/^'(e) ^ 00 as e ^ 0. In partic- 
ular, if we take /i^ = ^M^e~'^ , Lemma 11.31 implies that the domain {Re(z) < /x^} meets the 
pseudospectrum of iJ^ as e — > in the sense of Definition II. 2i It is also interesting to test the 
conclusions of Lemma ll. II in this situation: we know from ii) that ||e~*^'|| < C e"'^'* for some 
C > 1, but it follows from i) that C > exp(//e/^'(e)). Thus, although we can prove that 
the semigroup e~*^'= decays with an exponential rate iJ.^ ^ ^(e), this estimate is inadequate for 
small times because the prefactor C is exponentially large. In contrast, if we choose /i^ < K^'(e) 
for some k < 1, we can take C = ©(e^^) by Lemma O-iii). 

Proposition 11.91 will be proved in Section [5] below, using a complex deformation method 
(which exploits the analyticity properties of /) and the same localization techniques as in the 
proof of Theorem 1 1.81 We shall also give heuristical arguments, supported by accurate numerical 
computations, which show that the lower bound in (11.170 is optimal in the sense that the 
exponent v' cannot be improved. 

The rest of this paper is organized as follows. In Section [2] we prove some general properties 
of ^'(e) using mainly compactness arguments. Section [3] is devoted to the proof of Theorem 11.51 
using the approach developed by Villani [23]. We also show how that the variational method 
can be adapted to give the lower bound on ^'(e) in Theorem 11.81 up to a logarithmic correction. 
In Section m we give accurate bounds on the resolvent (H^ — iX)^^ for A G M using a dyadic 
decomposition of the real axis and semiclassical subelliptic estimates. In particular we prove 
Theorem 11.81 and we also provide a family (ci;(;)^g(o,i] of complex domains which avoid the 
pseudospectrum of as e — > according to Definition 11.21 with some non trivial e-dependent 
geometry. Finally, in Section [5l the lower bound (jl.l7p on S(e) is proved and illustrated by 
numerical computations. The proofs of Lemmas II. H [L3l and 11.71 are collected in Appendix lAl 

Acknowledgements. The authors would like to thank V. Bach, Y. Colin de Verdiere, M. 
Hitrik, A. Joye, G. Perelman, C. Villani, S. Vu Ngoc, and M. Zworski for fruitful discussions. 



2 Compactness estimates 

In this section we establish a few general properties of the quantities S(e) and ^'(e) defined 
in (|1.7p . (jl.Sp . We first show that any nontrivial function / does affect the spectrum of the 
operator H^. 

Lemma 2.1 Assume that f G L°°(M) is not a constant. Then 11(e) > ^'(e) > 1 for all e > 0. 

Proof: Suppose on the contrary that ^(e) = 1 for some e > 0. Then, in view of ()1.8p . (II. 9p . 
there exists Aq G M such that IK-ffe — zAo)~^|| = 1. Thus, for any n E N*, we can find G L^(M) 
with \\(j)n\\ = 1 such that 

1- - < \\{H,-iXoy^^n\\ < 1 . 
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Let ipn = {H^ — iXo) ^(pn, so that ipn & 'C' and cpn = (H^ — iXo)ipn- For all n G N* we have 

By Rellich's theorem, the sequence {ipn} is compact in L^(M). Thus, after extracting a subse- 
quence, we can assume that ipn converges in L^(IR) to some limit ip. By construction \\ip\\ = 1 
and 

\\ip'f + \\xibf < liminfdKf + llxV-nf ) < 1 . 

n— >oo 

Since + ll^J'i/'lP = {Hooip,ip), this inequality implies that ip{x) = vr"-^/^ e~^'^/^ (up to a 

constant phase factor). Moreover, since ipn^ipin L2(R) and \\{H^ - iXo)ipn\\ < 1 for ah n, we 
necessarily have \\{H^ — iAo)V'|| < 1- But {H^ — iXo)ip = ip + i{€~^f — Xo)ip and therefore 

\\{H,-zXo)iPf = y^|^(x)|2{l+(-/(x)-Ao) ]dx = l + y^|V'(x)|2(-/(x)-Ao) dx , 

hence it is possible to have — iAo)V'|| < 1 only if /(x) = eAo almost everywhere. □ 

Using similar arguments, we next prove that ^'(e) ^ oo as e ^ if and only if the level sets 
of / have empty interior (we assume here, for simplicity, that / is continuous). 

Proof of Proposition II. 4t Assume first that ^'(e)~^ does not converge to zero as e ^ 0. 
Then, according to definition (jl.Sp . there exist a positive constant 5, a sequence {e„} of positive 
numbers, a sequence {A„} of real numbers, and a sequence of normalized vectors in L^(M) 
such that e„ — > as ra — > oo and 

||(i?,„ - iXny^UnW > 6 > , for all n G N . 

Let Vn = {H^^ — iXn)~^Un, so that S P and Un = {H^^ — iXn)vn- The sequence {vn} is 
bounded in L^(M), because \\vn\\ < ||^n|| 1 for all n G N. Moreover, we have the identity 

{Un,Vn) = Wv'nW'^ + WxynW"^ + — if Vn,Vn) - iXnWVnW'^ . (2.1) 

Taking the real parts of both sides, we find + ||xt;„|p < \{un,Vn)\ < 1 for all n G N. 

By Rellich's theorem, the sequence {vn} is therefore compact in L^(M). Thus, after extracting 
a subsequence, we can assume that Vn converges in L^(]R) to some limit v. By construction, 
V G ^^(R) and \\v\\l2 > 6 > 0. 

On the other hand, multiplying both sides of (|2.ip by e„ and taking imaginary parts, we see 
that 

{fvn,Vn) -enA„||w„|p = enlm(n„,f„) > . 

Since {fvn,Vn) — > {fv,v), we conclude that 

{fv,v) def 

hm e„A„ = .. ..^ = fi ■ 

n—>oo \\V\\ 

Finally, for any test function (p G C^(M), we have 

en{Un,(p) = en{{H^„ - iXn)Vn,(p) 

= en{Vn,Hoo(p) +i{fVn,(p) - ie„A„(t'„, (/)) . 
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Taking the limit n — > oo, we obtain {fv,(f)) = fj,{v,(f)) for all (f) G (7q°(M), which implies that 
(/ — n)v = 0. Since / and v are continuous functions and v ^ 0, this is possible only if f~^{fJ.) 
has non-empty interior. Thus we have shown that, if all level sets of / have empty interior, then 
necessarily ^'(e) ^ oo as e ^ 0. 

Conversely, assume that there exists ^ G M such that f~^{fi) contains a nonempty open 
interval /. Choose cp £ C^{I) such that </> ^ 0, and let 

^ = HoocI) = {He - (for any e > 0) . 

Observe that ^ 0, because ^ a{Hoo). Thus 

^e)-' > IK^e-^r^ll > 
hence ^'(e) < ||V'||/||<^|| for all e > 0. □ 



3 Variational estimates 

The purpose of this section is to adapt the general method developed in [24] to our concrete 
problem, and to show that it allows to obtain precise estimates on the semigroup and the 
resolvent of the operator H^. To demonstrate the efficiency of this approach, we first give a short 
proof of Theorem 11.51 in Section 13.11 We next consider in more detail three model problems, 
and we show how to modify the proof of Theorem 11.51 so as to obtain optimal estimates in each 
case. With this information at hand, we consider in Section [3.31 the general situation where the 
function / satisfies Hypothesis [L6] and we prove the first half of Theorem [TTHl namely the lower 
bound on ^'(e) up to a logarithmic correction. We mention however that a complete proof of 
Theorem 11.81 will be given in Section HI using a different method. 

3.1 Proof of Theorem [TSl 

This section is devoted to the proof of Theorem ll.5i Consider any / G C^(M) such that /" and 
belong to L°°(M), and assume that estimate (II. lip holds for some v G (0, 1/2]. Let u{x,t) be a 
solution to the parabolic equation (jl.2p . namely 

dtu{x,t) = d'^u{x,t) — x'^u{x,t) f{x)u{x,t) . (3-1) 

To control the evolution of u{x,t), we introduce the quadratic functional 

m = ^(i|«|2 + |(|5,n|2 + x>|2) + /?Re((a,i2)i/'(x)u) + ^/'(x)>|2)dx, (3.2) 
where a, (3, 7 are positive constants to be determined below. We assume that 4/?2 < 07, so that 
m < I^Qul' + ^{\dM' + x'\u\') + ^f'{xf\u\')dx . (3.3) 

Note that, according to the general method introduced in [24J, the functional <I> is a linear 
combination of the quantities ||n|p, ||An|p, Ke{Au,Cu) and ||Cii|p, where as in the introduc- 
tion A = dx + X and C = [A,B] with B = {i/e)f{x). To compute the time derivative of <I>(t), 
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we use the identities: 



/ |updx = Re / udtudx = — {\dxu\'^ + x'^\u\'^) dx , (3.4) 
Jr Jr Jr 

/ a(|f?i.up + dx = — a|c?^ti — x^up dx — Re / a{dxu)-f'{x)udx , (3.5) 

Jr Jr Jr e 



1 d_ 
2di 
1 d_ 
2dt 

^ Re /" /3(9^n)i/'(x)udx = -- [ (3f{xf\u\'^dx-Re [ I3uif"' {x)dxudx 

+ 2Re / /3(a^n)i/'(x)(a^n-x^n)dx , (3.6) 

JR 

~ I 7/'(x)Vl'dx = - / 7/'(x)2(|5.^z|2+xVl')dx 
^ Clr Jffi Jr 

-2Re / 7/'(x)/"(x)tic>,ndx, (3.7) 
Jr 

which follow easily from (j3.ip after some integrations by parts. To estimate the various terms 
in the expression of ^'{t), we denote 

Kj = sup|(9^/(x)| , j = 2,3 , 

and we use the following bounds: 

-aRe I {dxu)-f{x)udx <- [ \dxu\^ dx + ^ [ f'{xf\u\^dx, (3.8) 
Jr e 4 7jj 7r 

-PRe [ uif"'{x)dxudx < K3P [ \udxu\dx < K3P [ (la^nj^ + x^l^) dx , (3.9) 
Jr Jr Jr 

2/3Re [ idxu)if\x){dlu-x\)dx < ^ [ \dlu-x\\'^ dx + — [ f'(xf\dxu\'^ dx, (3.10) 
Jr 2 Jig a J^ 

-27Re [ f'{x)f"{x)udxudx < \ f jO^up dx + 47^^! f /'(x)2|npdx. (3.11) 
Jr 4 Jig Jk 

We thus obtain 

^'{t) < (-1 + 1 + 1 + ^3/3) / (|a,u|2 + xVl')dx 
^ 4 4 /Jig 

+ (-a + I) ^ 1^2^ - x\|2 dx + (-^ + ^ + 47^K|) ^ nx f\u\' dx 
+ (-7 + ^) / /'(x)2(|5,n|2 + x>|2)dx . (3.12) 



We now fix 



Then 4/?^ < a7, 47^ ETl < /?/(4e), and we arrive at the simpler estimate 

<^'{t) < -1 / {\dxu\^ + x^\u\^)dx- ^ [ f'{xf\u\'^dx (3.14) 
4 Jr 2e Jig 

■ / la^u-x^pdx- J / /'(x)2(|a^M|2 + x2|n|^)dx . (3.15) 
Jr 2 Jig 
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In what follows we neglect both terms in (j3.15p and keep only the upper bound (j3.14p . Combining 
()3.14p with (jl.lip and assuming (without loss of generality) that e < 2/3, we arrive at 

^'{t) < -^(|a.np + x2H2)dx-^^/'(x)2H2dx-^(^)'^^H2dx. 

Using finally (jS.Sp . we conclude that < — r/<I>(t), where 



1 P Ml /2P 



Since < 1/2 and a,/3,7 are given by (|3.13p . there exists M2 > such that tj > 2M2€~'^ for all 
eG (0,1]. 

What we have done so far is to construct a quadratic functional ^ which is equivalent (with 
e-dependent constants) to the square of the norm in 

and which satisfies $(t) < e~''(*~*^$(s) for all i > s > 0. To deduce some information on the 
semigroup e~*^'= in L^(IR), we use the following standard argument. Let u{x,t) be a solution 
to p.ip with initial data uq G L^(R). In view of (I3.4p . we have 

^\\\dMt)f + \\Mt)f)dt = i(|Kf -IKOf) < ^iKf , 

hence there exists r G (0, ^/e\ such that 

||n(r)||2 < llnof , ||9,n(r)f + ||xn(r)f < 



2^ 

Thus, using the definition (j3.2p of <I>, the values ()3.13p for a, (3, 7, and the fact that |/'(x)| < 
|/'(0)| + K2\x\ for all x G M, we see that 

C 

HV^) < < -hof , for some C > independent of e G (0, 1] . 

e 



\u 



{t + Ve)f < mt + V^) < 2e~''*$(Vi) < — e-^*||nof , 



As a consequence, for any t > 0, we obtain 

2C 

e 

hence 

||g-{t+v/i)/f. II < 2^Nj 1/2 ^ alH > . 

Since we also have ||e~*''^'=|| < 1 for all t > 0, it follows from Lemma ll.ll -i) that S(e) > M2/e'^. 
Moreover, proceeding as in the proof of that lemma (see Appendix |A]), we find 

^(e)-^ < y""||e-*^<||dt < dt + I'minjl, | 

Since < 1/2, we deduce that ^'(e)^^ < C'e'^ log(2/e). This concludes the proof of Theorem 1 1.51 

□ 
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3.2 Three model problems 

A simple and general result such as Theorem 11.51 cannot give optimal estimates for all possible 
choices of the function /. In particular, if / satisfies Hypothesis ll.6| the lower bounds (|1.12|) 
with u = {k + 2)~^ (see Lemma ll.7p are too pessimistic. In this paragraph we consider three 
representative examples, and we show in each case that the proof of Theorem 11.51 gives in fact an 
optimal lower bound on ^(e) provided we choose in an appropriate way the parameters a, P, 7 
entering in the definition (13. 2p of the functional <I>. These examples are the building blocks which 
will allow us to treat in Section [3.31 the general case of a function / satisfying Hypothesis 11.61 It 
is interesting to note that the proof of Theorem 11.81 in Section [4] involves a similar enumeration 
of cases, see Proposition 14.11 



3.2.1 The case of a nondegenerate critical point 

We first consider the situation where / has a unique, non nondegenerate critical point, and |/'| 
is bounded away from zero outside a neighborhood of this point. The simplest example of such 
a function is f{x) = . In that case we can just follow the proof of Theorem 11.51 and choose 
the parameters a, /?, 7 as in ()3.13p . Indeed, the lower bound (jl.ll|) holds with v = 1/2, hence 
it follows from (j3.16p that < —r}^(f) with 77 = ©(e^^/^). This result is clearly optimal, 

because we know that S(e) = Re((l + i/e)!/^) = OifT^I'^) if f {x) = x^. 



3.2.2 The case when / has no critical point 

We next study the case when \ f\ is bounded away from zero. The simplest paradigm here 
is fix) = X. In this example, we have 

= -dl^ x^ ^-x = -dl^ {x ^ ^ ^ , 

hence S(e) = 1 + (4e^)~^. One can also estimate the resolvent of i^e, and one finds that 
^(e) = 0{e-'^l'^), see Proposition ST) To recover this result using the variational method, it 
is necessary to modify slightly the proof of Theorem 11.51 Indeed, if we still choose a, /3, 7 as 
in (j3.13p . then ()3.16p again implies r] = 0(e~^/^), because the lower bound (jl.lip holds with 
1/ = 1. However, since /' = 1, it is more natural here to set 7 = in (13. 2p and to assume that 
< a/4, so that <&(t) is strictly positive. Then proceeding exactly as above and using the fact 
that the terms <\3.7h . (13. 9p . and ()3.1ip vanish identically, we see that (|3.12|) becomes 

< (-1 + 1 + —) [ i\dM''+xM')dx 
\ 4 a / 

+ (-a + l) I \dlu-x\\'dx + (-f + ^) / /'(x)Vl'dx. 
We need a^ < e/3/2, so we can choose for instance 

2 1 1 

a = ( — P , and = [ — ^ 
V16/ V32 

Then we easily obtain < —r]^{t) with 



2 , 



. / 1 2(3\ _ 
= mm I — , — 1 = 0{e s) . 

Remark: In the more general case when /' is not identically constant, but still bounded away 
from zero, we have to deal with the additional term which appears in the left-hand side of 
Since \f"'\ < C\f'\, this term is easily bounded as in p.Sp . with a/e replaced by /3. 
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3.2.3 The case of a critical point at infinity 



Finally, we consider the situation where / has no critical point, but f'{x) decreases algebraically 
as X — > ±00. We fix A; > and focus on the following example: 

fix) = r XGM, 

which nearly satisfies Hypothesis II. 6[ In particular, f'{x) = C'dxl^'^'^^) as |x| — > 00. The proof 
of Theorem 11.81 given in Section H] applies to this example, and shows that ^'(e) = 0{e~'^) where 
v = 2/{k + 4). On the other hand, the lower bound (jl.lip holds with v = [k + 2)~^ < D, hence 
the proof of Theorem 11.51 only gives ^'{t) < —ri^{t) with r] = 0{e~^). However, we shall show 
how to choose the parameters a, /?, 7 in (j3.2p to obtain the optimal result in this case too. 

So let us again consider the functional supposing as usual < 07/4 to ensure (|3.3p . 
and let us go through estimates (j3.8p to (j3.1ip . The main advantage we have here compared to 
the general situation considered in Theorem 11.51 is that \f"\ + \f"'\ < Cfc|/'|, where denotes 
(here and below) a generic constant depending only on k. Estimates ()3.8p and (j3.10p . which 
involve only /', are unchanged. Estimate (13. 9p is replaced by 

-/3Re / uif"{x)d^uAx <] [ \d^u\'^ dx + CkP'^ [ f'ixf\u\'^dx. 



4 

Finally, after integrating by parts, we bound the left-hand side of (|3.1ip as follows: 

-27 Re / f'{x)f"{x)ud,udx = 7 / {f"{xf + f'{x)f"'{x))\u\^dx < Ckl I dx . 

JR Jr. Jr 

Summarizing, estimate (|3.12p becomes 

^'{t) <-\f (|9,.n|2+^2|^|2)dx+ (-^ + 4+Cfc(/?2+7)) / f{x)M^dx 



/ \o^u — x u\ dx + ( —'y -\ / 

Jr ^ a ^ Jm 



/'(x)2(|a^n|2 + xVP)dx . 



We now assume that < /?e/4 and Cfc(/3^ + 7) < /3/(4e). Using in addition the lower 
bound (jl.lip . which holds for the value v = (A; + 2)~^ in view of Lemma ll.7^ we obtain 

<^'{t) <--!' (|5,n|2 + xVP)dx-^ / /'(x)>|2dx-^(^)'' / \u\''dx. 
4 Jr 4e J]j 4 V e / Jk 

It follows that ^'{t) < -?7$(t) with 



T] = mm 



3a ' 3e7 ' 2 V e 



1 2 k 2k+2 2 

If we now choose a = 2^''+'^ , /? = e '=+■1 , and 7 = 8e , we obtain r] = 0{e ) which is the 
desired result. 

3.3 Improved decay rate under Hypothesis 11.61 

Let us now consider an arbitrary function / satisfying Hvpothesis ll.6l The analysis of Section 
suggests considering three different regions of space (near the critical points of /, near infinity. 
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or away from those regions), and shows how to choose the parameters a, f3, 7 so that the 
functional ^ defined by (13. 2p decays with the optimal rate in each region. Since different choices 
are needed in different regions, it is very natural to generalize the definition ()3.2p and to allow 
the parameters a, /3, 7 to depend on the space variable x. So we consider again the functional 

m = + ^i\d.n\^ + ^VP) + /?Re((9.n) if'{x)u) + lf'{xf\uf) dx , 

where now a, (3, 7 are positive functions of x, depending also on e. To choose these functions 
in the simplest possible way, we first decide that 

a = , and 7 = ^ ( y) ' ' ^^.l?) 

Indeed this relation between 0,7 and /?, e was already used in (j3.13p . so if we impose (|3.17p our 
set of parameters will be convenient near the critical points of /, provided /? is small enough. 
Moreover, as is easily verified, the parameters we chose in Section 13.2.31 also satisfy (|3.17p . 
so this definition will be adapted near infinity provided P is of the order of e~''^^^~^^^ in that 
region. Finally, in the intermediate region we can choose the same parameters a, /?, 7 as near 
the critical points of /, because (as was observed in Section r3.2.2p the functional $ will decay 
with the suboptimal rate r] = 0(e~^/^) which is still better than what we have near infinity. 
Summarizing, the choice (j3.17p is appropriate in all regions provided /? is small near the critical 
points of / and (3 = C'(e~'^/('^+^)) near infinity. Note that ()3.17p implies 07 = 4/?^, so that the 
upper bound ()3.3p is still valid. 



(3{x) 



k 

' fc+4 



A 



Fig. 1: The graph of the function (3 defined in (pTTSj) . Here ^ Ae . 
that all critical points of / are contained in the interval [—A + 1, A — 1]. 



and j4 > is sufficiently large so 



It remains to choose the function (3 appropriately. Let /3o > be a small constant, and ^4 > 
be large enough so that that all critical points of / are contained in the interval [— A + 1, A — 1] . 
Both parameters are independent of e, and their precise values will be specified later. We define 
the function /3 : M ^ M+ by 



/3(x) 



/9o 
/3o 



A 



2k 



fe+4 



if \x\<A, 

if 1^1 < |x| < , 

if |x| > B^ , 



(3.18) 



14 



where /c > is as in Hypothesis 11.61 and = Ae 2{fc+4) _ 'pj^g graph of (3 is represented in Fig. 1. 
It is easily verified that /3 has the following useful properties: for any c > 0, there exists eo > 
such that, for any e G (0, €q) and all x G M, 

e^/3'{xf < c/3(x)i , e/3{x) < c , and e(3'{xf < c(3{x) . (3.19) 

We now compute the time derivative of taking into account the additional terms involv- 
ing derivatives of a, f3, 7 which result from integrations by parts. We do not give the full details 
here, but simply indicate the main modifications with respect to the corresponding calculations 
in Section 13.11 Equation (13. 4p is of course unchanged, so we start with (13. 5p . Since a now 
depends on x, the right-hand side of (|3.5p contains an additional term which can be bounded 
as follows: 

— Re / a'{dxu){d'^u — x'^u) dx < - f a\d'^u — x'^u\'^ dx + f — \dxu\'^ dx . 
Jr 4 Jjj Jr a 

Remark that, in view of (j3.19p . we have for e > small enough: 

\ ; = ; , , ,i < — , for all X G M . 

a{x) 8/5(x)3/2 - 12 ' 

On the other hand, since = /3e/4, we have by (|3.8p 

-Re / a{d^u)-f'{x)udx <- [ |9^npdx + — [ Pf'{x)'^\u\'^dx. 



We conclude that 

„2„,|2 



1 d 
2di 



xU — X u\ dx 



[ a{\dxu\^ + x^\u\')dx <- [ ^\dl 
Jr Jr 4 

+ 11 \d^u\^dx+ [ lf{xf\u\^dx. (3.20) 
3 Jr Jr 4e 

We next consider ()3.6p . In that case one should add to the right-hand side the terms 

Re / j3'uif' {x){d1u — x'^u) dx — Re / jS'uif" {x)dxudx . 
Jr Jr 

The first one is easily controlled, for e small enough, by 

Re / I3'uif{x){dlu-x^u)dx < / '^\dlu - x^u\^ dx + f '^f{xf\u\^ dx 
Jr Jr o Jr O! 

< [ ^\Qlu-x\\'dx+ [ ^f'{x)^\u\^dx , 
Jr o jr I2e 

because f3''^/a <C /3/e due to (j3.17p . (j3.19p . Similarly, the second term is estimated by 
-Re / fi'mf"{x)dxudx < ^ [ |9^updx+ / 3P'^f"{xf\u\'^dx 

< ^ [ \dxu\'dx+ [ JLf^^f^ul'dx, 
J-^ Jr Jr J-^e 
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where in the last inequality we have used (j3.19p . as well as the fact that \ f"(x)\ < Ck\f'{x)\ for 
all X in the support of /?'. On the other hand, we have as in ()3.10p : 

2Re / P{d^u)if'{x){dlu-x^u)dx < [ ^\dlu - x'^u\'^ dx + [ ^f'{xf\d^u\^dx, 
Jm. 2 J 2 

where we have used the fact that 4/3^ = aj. Finally, the term involving f"'{x) in (|3.5p is 
estimated as follows. We first note that 

-Re / puif'd^udx < ^ [ \d^u\'^dx + 3 [ /?V"VPdx + 3 / p"^ f"'^\u\^ dx . 
Jr 12 Jr J\x\<a J\x\>a 

In the region |x| < A, we have 

3 [ P^f"'{x)^\u\'^dx < Sp^Ki [ l-updx < SP^Kl [ (|a^n|2 + x^jup) dx , 
J\x\<A Jr Jr 

and we choose Pq > small enoug h so that 3(3^Kl < 1/6. In the region |x| > A we know that 
\f"'(x)\ < Ck\f'{x)\, and in view of (j3.19p we find for e small enough: 

3 / P^f"'ixf\u\^dx < [ J^f'{xf\u\^dx. 
J\x\>A Jr 12e 

Summarizing, we have shown that 

Re / p{d^u)if'{x)udx < - [ ^ f'{xf\ufdx + l [ + x>|2) dx 

dt Jm 4e 3 

+ / ^\dlu - x\\^ dx + [ ^f'{xf\d^u\^dx . (3.21) 
Jr ° 2 

Finally we turn our attention to (j3.7p . This time the right-hand side contains only one 
additional term, which can be estimated as follows: 

-Re / -r'f'ixfuda^udx < [ ^f'{xf\d^ufdx+ [ ^f'{xf\u\^dx 
Jr Jr 4 7 

< / lf\xf\dM'dx+ [ ^f\xf\u\^dx, 

where we have used the fact that 7'^/7 <^ (3/e if e is small, see ()3.17p . ()3.19p . It remains to 
bound the last term in the right-hand side of (|3.7p . Using a standard (e- independent) partition 
of unity and we decompose 7 as 7 = 71 + 72, where 71 and 72 are nonnegative Lipschitz functions 
satisfying supp(7i) C {|x| < 2A} and supp(72) C {|x| > 3yl/2}. Then we observe that 

-2 Re / 7i/'(x)r(x)na,txdx < / ^f'{xf\u\^dx + l2Kl [ 'S^ld^u]^ dx 
Jr Jr J\x\<2A p 

< [ r^nxf\u\^dx+^ [ \d,u\'dx, 

where in the second inequality we have used the relation = 64/3^ and the fact that 

P{x) < f3Q2^^ if |x| < 2A. We have also chosen /3o > small enough so that 768 K| 2^^ < 1/12. 
On the other hand, integrating by parts, we find 

-2Re / ^2nx)nx)udxudx = [ ^,{f"^ + f f")\uf dx + [ ij'f"\u\'dx . 
Jr Jr Jr 
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We know that \ f"\ + \f"'\ < Ck\f'\ on the support of 72. Moreover, since I7V7I < ||/37/3| < C 
by (|3.18|) . it is clear that I72I < C7 for some C > independent of e. Thus 

-2 Re / -f2f'{x)f"{x)ud^udx < C I -ff\xf\u\'^dx < I -^/'(x)Vl^dx , 

because 7 ^ /?/e if e is small, see (|3.17p . (I3.19p . We conclude that 

lA [ jf'^\u\'dx < - [ ^f{xf{\d,u\^ + x^\u\')dx 

+ ^ [ \dM^dx+ [ ^f{xf\u\^dx. (3.22) 
J-^ jR JR 4e 

Summarizing, if we collect estimates ()3.4p , ()3.20p , (j3.2ip , and (j3.22p , we obtain 

^'(t) < - 7 / {\d^uf + x'^\u\^) dx - [ j-f'{xf\ufdx (3.23) 
4 Jr Jr 4e 

- / ^\dlu - x\\^ dx - [ y'ixf{\d^u\^ + x^\u\^)dx , 
JR 8 JR 4 

and as in Section [3. II it is sufficient to keep only the first line in (I3.23[) . On the other hand, with 
our choice of the function it is proved in Lemma lA. II that there exists Mq > such that 

[ + x^\u\^ + ^ f'{xf\u\A dx>^ \\u\\l, , Where u = • (3.24) 

JR ^ e y e'^ A; + 4 

Combining ()3.24p with the first line of (j3.23p . we easily obtain 

cl,'(t) < -- f {\d^,uf + x^\u\^)dx- f llf'(xf\u\'dx-^ [ \u\^dx. 

Using now the upper bound (|3.3p . we deduce that ^'{t) < —r]^{t) if 

1 1 



r] = mm 



6||a||_Loo 6e 



As in Section \37l] we conclude that ^'(e)~^ < Ce'^log(2/e) for some C > independent of e. 
Up to a logarithmic correction, this proves the lower bound in (jl.l5p with the rate e^, which is 
better than e'^ since v = {k + 2)~^ < i?. 



4 Resolvent estimates 

In this section we obtain precise estimates on the resolvent of along the imaginary axis, 
and thereby prove Theorem 11.81 In doing so, we construct a nontrivial complex domain uj^, 
parametrized by e G (0,1], which avoids the pseudospectrum of as e ^ in the sense 
of Definition 11.21 The boundary of this domain is a graph over the imaginary axis which 
exhibits a complicated structure involving various e-dependent scales, see Fig. 2 below. This 
phenomenon does not occur in familiar examples such as self-adjoint operators, kinetic Fokker- 
Planck operators, or complex harmonic oscillators, see [4]. Our proof will show that the value of 
the quantity ^'(e) results from the competition between various microlocal models related either 
to the critical points of / or to the behavior of / at infinity. 
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Throughout this section, we assume that / : R ^ M satisfies Hypothesis 11.61 In particular 



/(M) = /(M) U {0} and / has only a finite number of critical points. The finite set of critical 
values of / is denoted by 

cv(/) = {/(x); xGM, /'(x) =0} . 

For any A G M and any e G (0, 1), we define 

K{e,X) = \\{H,-iX)-^\\ . (4.1) 
The following proposition gives accurate bounds on ^(e, A) in various parameter regimes: 

Proposition 4.1 For e £ (0,1) and A G the quantity ^(e. A) defined in ()4.ip satisfies the 
following estimates: 

i) //dist(eA,/(M)) >6>0, then ^(e. A) < e/6. 

ii) //dist(eA,cv(/) U {0}) > 6 >0, then ^(e. A) < C^e^/^. 

ill) If X = A(e) is such that lime^o ^^(c) = c G cv(/) \ {0}, then limsup^_^o e^"'^/^K(e, A(e)) < C. 
iv) For A = 0, the quantity n{e, 0) satisfies 



k{€,0) < < 



CeT^ if 0^ 

^f OG/(M)\cv(/) , 
_ Ce'"'"^^'^} if Oe cv(/) . 



v) There exists C > I such that n{e, X) < Ce''+4 for all (e. A) G (0,1) x M. Moreover, if 
'^(6) A) > C~ e'^+^j t/ien A zs comparable to e ^=+4 . 



2 2 

^'(e) = minK(e,A)"^ G M^^e"'^ , M4e~*+4 



Finally all estimates in i), ii), Hi), iv), and v) are optimal, in the sense that one can find 
X = A(e) so that the pair (e, A(e)) satisfies the required conditions as e — > and so that K{e, A(e)) 
is comparable to the upper bound in this limit. 

Theorem 11.81 is of course a direct consequence of Proposition 14. 1^ since 

for some M4 > 1 . 

Proposition 14.11 also allows to localize the pseudospectrum of accurately: 
Corollary 4.2 The complex domain uj^ defined for e G (0, 1] by 

CO, = [j^zeC;\z-iX\< ^J' I U |z G C ; Re(z) < o| 
AeM ' 
avoids the pseudospectrum of as e ^ 0, whereas the domain 

We = G C ; Re(z) < ^'(e) \og{e'^f] n (C \ 

meets the pseudospectrum of as e ^ 0. 

Proof: If \z — iX\ < (2K;(e, A))^-*^, the resolvent formula shows that \\{H^ — z)"^ || < 2K(e, A) < C. 
Moreover, we know that \\{H^ — z)~^ || < 1 if Re(z) < 0, thus to^ avoids the pseudospectrum of 
as e — > 0. Remark that u;^ contains the half-plane {z £ C; Re(2:) < ^{e)/2}. On the other hand, 
if = ^(e)log(e~^)2, it follows from Proposition g?!] that > ^'(e)(l + log^'(e) + log(e~^)) 
if e is small, hence the half-plane {Re(z) < /i^} meets the pseudospectrum of by Lemma ll.3l 
The same conclusion holds of course for ojf. □ 
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Fig. 2: The domain uj^ on the left-hand side of the solid curve avoids the pseudospectrum of as e —* 0. The 
picture on the right shows the geometry at the scale ez = 0(1), while the left picture focuses on the region where 
ez is small. Here = {2 € C ; Rez > 0, min/ < elm 2; < max/} and cv{f) — {ci, 02,03}. 



To prove Proposition I4.H we start with a variant of the so-called IMS-localization formula, 
see [2J. 

Lemma 4.3 Consider a Schrddinger-type operator Q = —A + V in M.'^, where V is a measurable 
function. Take a locally finite partition of unity {xj}jej, where Xj G C^(IR'^,M), such that 

^Xj{xf = l, forallxGR'', (4.2) 

and 

ml = sup ^|Vxj(x)p < +00 , ml = sup ^(Axj(x))^ < +00 . (4.3) 
Then the estimate 

2\\Qu\\'^ + 3ml\\uf + 8ml\\Vuf > ^ HQxj^^f , (4.4) 

holds for any u G Cq^{W^). Moreover, ifYieV > 0, then 

2\\Quf + Zml\\uf + SmlYle{Qu,u) > "^WQxjuf ■ (4.5) 

Proof: For any x e C^(M'^,R) we have 

Q*x'Q = Q*xQx + Q*x[x,Q] 

= xQ*Qx + [Q*,x]Qx + Q*x[x,Q] 

= xQ*Qx + [Q* , x\xQ + Q*x[x,Q] + [Q* ,x][Q, x] 

= xQ*Qx - [A, x]xQ + Q*x[A, x] + [A, x] [A, x] ■ 
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Since 

[A,x] = 2(Vx) • V + (Ax) = 2V • (Vx) - (Ax) , 

this implies 

Q*x^Q = xQ*Qx - V • {Vx^)Q + (Ax)xQ + Q*(Vx') • v + Q*x(Ax) - KR 



with = 2(Vx) • V + (Ax). We now apply this identity with x = Xj and sum over j £ J. 
In view of (j4.2p . the left-hand side reduces to Q*Q, and the second and fourth terms in the 
right-hand side disappear, so that 



Q*Q = Y.XjQ*Qxj + Y.[(^^^)xjQ + Q*x^(^x^))-Y.Vx, ■ 

jeJ j€J jeJ 

Thus, for any u e C^(M'^,M), we have 

WQuf > ^{WQXjuf - iWxjQuf + ||(Ax,)nf ) - (8||(Vxi) • Vtxf + 2\\iAxM\')) 



which is (j4.4p . When ReV >0, the inequality ||Vu|| < Re{Qu,u) implies ()4.5p . □ 

The idea is now to apply Lemma 14.31 to the one-dimensional operator H^: — i\, using a 
dyadic partition of unity. This allows us to reduce a global problem on the whole real line to 
a family of compactly supported problems, indexed by the parameter j G N. The choice of a 
dyadic partition is convenient to take into account the precise behavior of the coefficients of 
as X — > ±00. 

Lemma 4.4 For j G N, e > 0, and A G M, consider the operator 

P,- = + 22^x2 + '-f{2^x) - iX , (4.6) 

and let 



Cj{e,X) = inf|||Pj-,,A^i|| lueC^iR), suppn C Kj , ||m|| = l} , (4.7) 

where Kq = [—1, 1] and Kj = [—1, —1/4] U [1/4, 1] for any j > 0. Then the quantity K{e, A) = 
\\{H^ — iX)~^\\ satisfies 

inf C7,(e,A)V^ < «(e,A) < cf inf C7,(e, A) , (4-8) 
for some constant C > 1 independent of€,X. 

Remark 4.5 It is clear that Cj{e, A) > 1 for all j G N, e > 0, A G M, because 

< Re {Pj^^^\U ,u) < ||Pj.e,A^|| ll""!! ) 

for allue C^{R). 
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Proof: We choose a dyadic partition of unity {Xj}jeN such that 

oo oo „ 

j=0 j=l 

where XOjX ^ C'o°(M) satisfy 

, . / 1 if kl < I , , , , / 1 if i < |x| < I , 
' I if |x| > 1 , ' I if |x| < I or |x| > 1 . 

It is clear that such a partition exists and that the quantities ml, defined by ()4.3p are finite. 
Thus we can apply Lemma 14.31 to the operator Q = — iX = —d"^ + + ~f{^) ~ ^-^i for 
any e > and A G M. Since \\uf < Re{Qu,u) < \\Qu\\ \\u\\ < \\Quf for all u G C;f (M), it 
follows from (14.51) that 



j=0 

where = 2 + Smf + Smg. Now, for any j G N, we define 

Vj{x) = 2^/^Xji'^^x)u{2^x) , xGR , 

so that suppvj C suppxj(2-'-) C Kj and {Pj^^\Vj){x) = 2^^'^{Qxju){2^ x). If we denote m(e, A) = 
infjgpj Cj{e, A), we thus find 

oo oo 

C^WQuf > ^ ||P,-e,At^if > m{e,Xf^\\vjf > m(e,A)^||nf , (4.9) 

j=0 j=0 

because J2j IbilP = J2j IIXj^lP = II'^^IP- Since Q = — iX and Cq°(M) is dense in L^(]R), it 
follows from (gJD that ^(e, A) = \\{H^ - iX)-'^\\ < C/m{e, A), which is the upper bound in (j4.8p . 

To prove the lower bound, we fix 5 > 0, e > 0, A G M, and we take j G N, Vj G C^(M) such 
that Vj ^ 0, suppiJj C Kj, and ||Pj^e^A^^j|| < ("^-(e. A) + Setting Vj{x) = 2^ /'^u{2^ x), we 

find that \\Qu\\ < (m(e, A) + 5)||tt||, hence K(e, A) > (m(e. A) + 6)~^. Since (5 > was arbitrary, 
this concludes the proof. □ 

Proof of Proposition I4.lt 

The proof of assertion i) is easy: if dist(eA, /(M)) > 5, then dist(iA, 0(ife)) > S/^, where @{H^) 
is the numerical range defined in (jl.Sp . The last inequality in (jl.9p then implies K{e,X) < e/6. 

In the remaining four cases, we start from Lemma 14.41 and use (|4.8p to bound K{e,X). It 
turns out to be convenient to rewrite the operator Pj^e,x in the equivalent form 



.e2('=-2)j^2 ^ g2('=+2)ix2 + i(2^J/(2^x) - 2'=J'eA)l , (4.10) 



where A: > is the parameter that governs the asymptotic behavior of f{x) as |x| — > oc, see 
Hypothesis 11.61 For later use, we observe that 

Cf =^ sup sup 2''^\f{2^x)\ < +00 . (4.11) 
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ii) Suppose that dist(eA, cv(/) U {0}) > 5. Without loss of generahty, we also assume that 
< because otherwise we can use the estimate established in i). For any u £ C^(M) 

with suppu € Kj and u^O, we have the lower bound 

\\Pj,e,xu\\ |Im(Pj-,.An,n)| ^ |([2^i/(2J-.)-2feJ-eA]n,n)| 1/ . . _ 

\\u\\ - ||n||2 e2fci||'u||2 " e ' 2^^ J ' 

where C/ is given by (j4.1ip . Since e|A| > 6, we deduce that Cj{e,X) > 5/{2e) whenever j is 
large enough so that 2^^ > 2Cf /6. Thus only a finite number of indices j have to be considered, 
and the problem is therefore reduced to finding a lower bound on the quantity \\{H^ — iA)M|| 
when u G C^{{x G R; |x| < Rs}), for some Rs > 0. On a bounded domain, we can drop the 
bounded term in and only consider the operator Q = —d'^ + f (/(x) — eA). Take two 
functions 6*0, Oi £ C°°(M,]R) such that supp6'o C [-2,2], supp6'i C (-cx),-l] U [l,+oo), and 
Oq{x)'^ + 9i{x)^ = 1 for all x G M. Given a > which will be fixed below, we consider the new 
partition of unity 

Xoixf + Xi{xf = I , with Xj{oo) = for j = 0, 1 . (4.12) 



The quantities ()4.3I) for this partition satisfy mf < Cie and rrig < Cie ^'^ for some Ci > 0, 
hence ()4.5p implies 

2\\Quf + ^\\uf + -^WQuW M > \\Q{xou)f + \\Q{xiu)f . (4.13) 
The last term is easily bounded from below using the support condition on xi, which yields 

||Xi«IIIIQ(Xi^)ll > -m{x)-\€)xiu,xiu)\ > e'-^Xi^f • (4.14) 
e 

To bound Q{xou), we observe that f^^{eX) is a finite set {xi, . . . ,x„} because / has a only a 
finite number of critical points. Since dist(eA, cv(/) U {0}) > 5, it follows that for e > small 
enough the support of xo is the disjoint union of n intervals /i, . . . , /„ of size 0{e") centered at 
the points xi, . . . , x„. In particular, xo"" = l^"=i '^t with suppn^ n suppn^/ = when i ^ i' . 
Inside I^, the operator Q is well approximated by = —d^ + -f'{xi){x — xg), because 



f{x) - f{x() - f'{xi){x - Xi) 



e 



> JiiQ.-.f-^^^^l^ikir (4.15) 

On the other hand, the operator Qi is unitarily equivalent to the microlocal model 

= jy^{-d^±iy) , where 7 = , (4.16) 



57 
since 



\Puf = \\d^uf + \\yuf ± 2i{u , dyu) > -\\dyuf + \\yuf - C\\u\\ 



which satisfies ||(57n|| > C72/^||ii||. Actually, the operator P = —dy±iy has a compact resolvent, 

with an empty kernel {Pu = implies = Ke{Pu , u) = 0), and it is therefore invertible. 

We refer to Chap. 27 in [16] or to [4] for a geometric analysis of similar and more general 
models. Alternatively, one can use the Lie-algebra approach developed in [20], [Tlj if one sets 
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—dy + ydt = —Xl + Xq with Xi = dy, Xq = ydt and if one considers the value r = 1 for the 
frequency variable r dual to t G M. Summarizing, we have shown that 



\\Q{xou)f = ^WQuef > [Cse-^/^-C2\\f"\\l^e^'^-yxouf , (4.17) 
e=i 

for some C3 > 0. We now assume that a > 1/6, so that e^"""^ <C e~^/^ if e is small. Replacing 
(I4T4]) . ([4171) into (I4l3]) . we thus find 

C4 II ii9 . ^—1 . f 9o-— 9 —A/^^ 'lulP 



C4||Qu|p + ^||uf > C4-imin{62--2,e-V3| 



for some C4 > 0. Finally, we suppose that a < 1/3, so that e ^ e ^ e^'^ ^. We thus 
obtain the estimate 2C|||Qm|P > e~^/^||u|p, which proves that ^(e. A) < Ce^/'^. 

iii) The assumption lim^^o = c G cv(/) \ {0} implies that e|A| > 5 for some fixed 5 > 
if e is small enough. Thus we can reduce the analysis to a bounded spatial domain like in 
case ii). By assumption /~^(c) is a finite set which contains at least one critical point of /, 
but in general this set contains noncritical points as well. Using a partition of unity, we can 
treat the noncritical points separately and estimate their contributions as in case ii). So, for 
simplicity, we assume here that /~^(c) consists of critical points only. We shall consider two 
different regimes, depending on how fast eA converges to c as e ^ 0, and then check after some 
iterations that they provide a complete information for the general assumption lim^^o e-^ = c. 
a) We first study the case when eA converges slowly to c. More precisely, we assume that 

e"' < \eX-c\ < (4.18) 

where < (T2 < < 1/2 and 3(72 > 5c7i — 1. If (T2 = 0, we also suppose that eA ^ c as e — > 0. 
Several intervals [(T2,o"i] will be fixed iteratively in step c) below. Our goal is to show that 
\\Qu\\ > Ce-^/^\\u\\ for ah u E C^{R), where Q = -dl + ^(/(x) - eA). We choose a G (ai, 1/2) 
such that 

2(7l 1 1 CJi (Ti — (72 , , 

and we use again the partition of unity defined by (j4.12|) . If x G supp xo oi x G supp x'l , we 
have |/(x) — eA| < 2e'^ and using (j4.18p we easily find 



^(Ti/2 I f(^\ _ „|l/2 

_ < 1-^^^^^ "I < |/'(x)| < Ci|/(x)-c|i/2 < C2e-/2 



C2 Ci 

for some Ci, C2 > 0. These estimates allow to give precise bounds on the quantities m\ and 
defined in (14.31). Since 



we obtain 



ml < Cse"^-^" , and mj < Cl{e-^'' + e^''^-*'') < 20^ e^"^-^'' , 
for some C3 > 0, and it follows from (j4.5p that 

3||Qnf + 22C|e2-^-4'^||uf > \\Qixon)f + \\Qixiu)f . (4.20) 



23 



By (|4.14p we have > e'^~^||xiM||, and to bound ||(5(xo^)|| we proceed as in case ii). 

Denoting f~^{eX) = {xi, . . . , we decompose as before xo"" = Z]"=i '^^ observe that 

every x G supp ug satisfies \x — xg I < 2C2e'^-'"i/^ because 

2e'^ > |/(x)-eA| > |x - inf{|/'(y)| ; y G supp^x^} > C^^ e'^^/^lx - xg\ . 

As in ()4.15p . we thus find 

for some C4 > 0, because e-^/^\f'{xE)\^/^ > C^'^^^e^'l > e^'^-2<7i-2 gjg])^ Summing 
over £ as in (j4.17p . we obtain the desired lower bound on ||(3(xo^)|Pi and returning to (j4.20p we 
arrive at 

for some C5 > 0. Since e2o-2-4cr ^ ^ ^2(7-2 (gjg])^ conclude that 

2Cl\\Quf > e^-t||nf > e~^||nf , for ah u G C^(M) . 

b) We now assume that |eA — c| < , for some a G (5) ^)) and we use again the partition of 
unity defined by (|4.12p . As before, we have < Ci\f{x) — c\^^'^ < C2 e*^/^ for all x G suppxo 

and all x G suppx'i, and it follows that the quantities ()4.3p satisfy ml < C3 e~'^, m| < C| e~^'^. 
Thus (j4.5p becomes 

3\\Quf + 22Cle-^^\\nf > \\Q{xou)f + \\Q{xiu)f > \\Q{xou)f + e^^~^\\xiuf. 

Since e"^'^ <C e"^ < e^'^"^, it is sufficient to show that ||Q(xo^^)|P > Ce'^WxouW^ ■ To this end, 
we consider the set /~^(c) = {xi, . . . ,Xn} and we decompose as before xou = Yl^=i the 
support of Ui, we can use the quadratic approximation 

Q, = -dl + 1 [^^^(x - x,f - (eA - c) 



e 
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Indeed, if x G suppti^, then \x — xi I < CAe"/^ for some C4 > 0, and it follows that 

for some C5 > 0. On the other hand, the operator Qi is unitarily equivalent to the microlocal 
model 

1/2/ a2 , ■ 2 ■ N u \f"{^i)\ J 1/2 eA-C 

^ ' {—o^. ±iy — i^jl) , where 7 = and /i7 ' = . 

^ 2e e 

Using the methods presented in [16j (Chap. 27), [llj, [4J, or even reproducing the analysis of ii) 
leading to (j4.16p . we find 

\\{-dl ± iy^ - i^JL)v\\ > C(l + |/i|^/^)||f || > C\\v\\ , for all v G C^(R) ■ 

Since e^'^^^ <^ e^^, this shows that HQn^p > Ce~-'-||ti£|p, hence also ||(5(xo'u)|P > C'e~"'"||xo^^|P- 



As in step a), we conclude that HQ^lP > C e for all u & C( 

c) Take = G (|, i) and, for any n G N, let cig = 0""+^ = ^af - i. Since ^ > 1 and since 
the solution to x = — | is x = | > cJ]', it is clear that (T2 < cr" for all n and (T2 —00 
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as n — > oo. Let no S N be the smallest integer for which aJ^ < 0. As ^ > |, the condition 
3o"2 > 5a" — 1 is satisfied for all n < hq. Applying step a) to all intervals [max{0, erg}, cr"] for 
n = 0, . . . , no, we obtain the lower bound \\{H^ — iX)u\\ > C e"^/^||'u|| whenever A satisfies ()4.18p 
with [(72, CTi] C [0, In other words, there exists a constant > such that 

e\{e)-c\>efo \ ^ ^imsup 6-V\(e, A(6)) < K 



lime^o eA(e) = c 



e-+0 



Meanwhile step b) provides the same conclusion (possibly with another constant K) under the 
assumption |eA(e) — c| < eso . This concludes the proof of case iii). 

iv) We consider the operator Pj.t^x given by (|4.1U|) for A = and i > 1. For any u G Cg°(M) 
with suppu C Kj = {\< \x\ < 1}, we have 

ll-ull ll-PjEOitll > \Re{Pjeou,u)\ > 2^^ / x^\u{x)f dx > 2^^~^\\uf , 
\\u\\\\Pj,,flu\\ > |Im(P,-,,on,^x)| > -^1^ 2^^\fi2^x)\\u\^dx > ^ \\uf , 



where mj = inf{2'^-^ |/(2-^x)| ; | < |x| < 1}. From Hypothesis 11.61 we know that rrij — > 1 as 
j ^ oo, so we can take J S N large enough so that mj > 1/2 whenever j > J. Thus 

C,(6, A) > 1(22^-4 + > C6-^ , for all j>J. 

The case when j G {0, . . . , J} corresponds to a bounded spatial domain, and can be treated 
exactly as in ii) and iii). We find that ||i^e^|| ^ C'e~'^||n|| where cr = 1 if /~^(0) = 0, a = 2/3 if 
E /(M) \ cv(/), and = 1/2 if G cv(/). 

2 

v) We want to show that K(e, A) < Ci 6^+4 for all (e. A) G (0, 1) x R. In view of Lemma we 
have to verify that 

Cj{e, A) > C2 , for all j G N and all (e. A) G (0, 1) x M . (4.21) 

As was already observed, the analysis of Cj{e, A) for j in a finite set {0, . . . , J} can be performed 
as in ii), iii) and yields (in the worst case) the lower bound Cj{e, A) > Cj e~^/^, which is already 
better than (j4.2ip . Thus is it sufficient to consider large values of j. Fix such a j and choose 
u G Co°(K) such that suppu C Kj = {i < |x| < 1}, \\u\\ = 1, and ||P;>,a'«II < 2Cj(e, A). 
Proceeding as in iv), we easily find 

> 2'^-' , and ||P,„.^|| > "^^^^^^"^^ , (4.22) 

where gj{x) = 2^^ f{2^x) — 2^^e\. If 2^ > e^'=+4, the first inequality already implies (j4.2ip . In 
the converse case, we integrate by parts and obtain the relation 

\\Pj,,,xuf + 2\\uf = ||Qj-,,Anf + 2||xu^f + 2^-''||x\||2 > \\Qj,,,xuf , (4.23) 

where Qj,e^x = Pj,e,\ — 2^^x'^ . Thus ||Pj^£_A^|| ^ ||Qj,e,A'w|| — V^ll""!!; and combining this estimate 
with (j4.22p we obtain the sharper result 

2C,(e,A) > ||P,-,,An|| > i (22^-4 + + ||Q^.^,^|| _ V2) . (4.24) 
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We next observe that 

Qi.^.A = ^ {-h^dl + ig,) , where = e2^^-^^^ = O(e^) . 

From Hypothesis 11.61 we know that the derivative g'j{x) = 2^^^^^^ f {2^ x) converges to —kx\x\~^ 
as j oo, uniformly on Kj. In particular, there exists C3 > such that > C3 and 

\gj{x)\ < C^^ for all X G Kj if j is sufficiently large. Thus, using a partition of unity and 
a Taylor expansion like in case ii), we can again reduce the analysis to the microfocal model 
(j4.16p . We obtain the following lower bound 

\\{-h^dl + igj)u\\ > C4 h'^^^\\u\\ , for all u G C^{R) with suppu C Kj . 

Returning to (j4.24p . we arrive at the estimate 

2C,ie,X) > + V2) >C,e-^^, (4.25) 

which proves (I4.2ip . 

2 

Assume now that there is a family (e. A) with A = A(e) such that e ^ and K{e, A) > Cq e'^+i 
for some Cq > 0. By Lemma 14.4^ there exists C7 > such that, for all (e. A), 

Cj{e, A) < C7 e"^ , for some j = j{e, A) G N . (4.26) 

It is clear that j{e, A) ^ cxd as e ^ 0, because if j stays in a finite set {0, . . . , J} we have 
the lower bound Cj{e,X) > Cje~^/'^ which is incompatible with (j4.26p . We also know that 
2^ < C e~'^^^^^^\ otherwise ()4.24p would contradict (j4.26p if e is small. Thus we can use estimate 
(j4.25p which implies, in view of (j4.26p . 

Cg^^e"^ < 2^ < Cge"^ , for some > 1 . (4.27) 
Next we choose u as before and we deduce from (j4.22p . (j4.26p that 

ini^eKj \9j{x)\ 



e. 



2kj 



2CTe~~i > 2Cj{e,X) > \\Pj,e,\u\\ > 
As gj{x) = 2^^ f{2^x) — 2^^eX, this implies by direct calculation 

inf 2^=^(2^'^) -2C7e-^ < A < ^ sup 2^=^(2^'^) + 2C7 e"^ . (4.28) 

Since 2^^f{2^x) — > uniformly in Kj as j — > c« by Hvpothesis 11.61 and since 2^ satisfies 

(fOT]) . it follows immediately from (f08]) that 

Cg e < A < Cg e ^+4 ^ 

for some Cg > 1, which is the desired estimate. 

Finally, it remains to check the optimality of the various estimates established so far. This is 
a rather straightforward task, because the arguments we gave to bound ^(e. A) from above also 
indicate how to choose an appropriate test function u G (M) which "saturates" the inequality 
and gives the corresponding lower bound. We shall provide the details only in the last case v), 
because this implies the upper bound on \I'(e) in Theorem 11.81 
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Fix xo S (1/4,1) and take j £ N very large. We define e > 0, A > 0, and /i > by the 
following relations 

2J = , h"^ = e2('=-2)j' , eA = f{2^xo) . (4.29) 

Next, we choose v £ C^(M) such that ||f || = 1 and suppi; C [—1, 1], and we denote 

Mx) = T^v(^^), xeR. (4.30) 



hl/3 V /i2/3 

It is clear that G Co°(M), \\ufi\\ = 1, and that suppu/j C Kj if /i > is sufficiently small, that 
is, if j is sufficiently large. Recalling that 

Pj,^,^ = -L{-h''dl + h'"^^'' + ig,{x)) , where g,{x) = 2^^f{2^x) - 2^^ eX , 



e2^i 

we claim that there exists Cq > independent of j (hence of e. A) such that 

\\Pi,e,xUh\\ <Co^ = Coe-^ . (4.31) 

This implies that Cj(e,A) < Coe'T^, hence K(e,A) > Cq'^ by K8\t . 

Checking (j4.3ip is of course straightforward. First, using (|4.30p . we find ||/i^c?^n/i|| = 
/j2/3||^//||_ jsjgxt, since x^ < XQ + 2|a:; — xo| for all x E Kj, we have \\x^Uh\\ < xl\\v\\ + 2h'^/^\\v\\ < C. 
Finally, since gj{xo) = by our choice of A, we have for all x G [1/4, 1] 

\gj{x)\ < \x-xo\ sup \g'j{z)\ < C|3;-xo| , 

^£[1/4,1] 

where the constant C does not depend on j by Hypothesis 11.61 Thus ||5jtt/i|| < C/i^/^||w||, and 
the proof of (I4.3ip is complete. □ 



K(e, A) 




Fig. 3: The resolvent K(e, A) = \\{H^ — iX) ^\\ is represented for e = 2 ^"^ and / given by (|4.32p . 
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To illustrate the various regimes appearing in Proposition 14. H we computed numerically the 
function ^(e, A) in the particular case when 

/(^) = (1 + ^2/3)3 ^ (4-32) 

Up to a multiplicative factor, this function satifies Hvpothesis ll.61 with = 4. It has exactly two 
critical values ci = 1, C2 = 27/16, which correspond to the critical points xi = and X2 = ±1 
respectively. The result for e = 2~^^ is displayed in Fig. 3. The large peak near the origin 
corresponds to the regime described in v) above, while the two sharp spikes are associated to 
the critical values of /, as explained in iii). The predictions of i), ii) are also confirmed, but 
the value e = 2~^^ is not small enough in this example to test the prediction of iv): we know 
from Proposition 14.11 that ^(e, 0) = 0(e^/^) since = 4 and ^ /(K), but in Fig. 3 the value 
computed for A = is still clearly below the ©(e^/^) spikes associated to the critical values of /. 



5 Spectral lower bounds 

This section is devoted to the detailed study of the particular example ()1.16p . We shall first prove 
Proposition 11.91 then illustrate the optimality of the lower bound (jl.lTp using numerical 
calculations and heuristic semiclassical arguments. As was mentioned in the introduction, the 
estimate (|1.17p implies that ||e~*-'^' || < C e"'^"* for some C > if /x^ = ^M^ , but it follows 
from Lemma [1.1 1 that C > e exp (/x,/^'(e)) where = Oie^-" ) ^ oo as e ^ 0. It is very 

likely that such an exponential bound in the semiclassical framework cannot be obtained while 
remaining in the C°° category. Thus, like in the study of resonances (see [T|[2|ll2tll5j). we shall 
exploit the analyticity properties of the function (jl.l6p by performing a complex deformation. In 
this example we shall observe an interesting phenomenon when comparing the pseudospectral 
and spectral estimates. In contrast to what happens for resonances, the large peak of the 
quantity ^(e. A) = \\{H^ — iX)~^\\ in the region A ~ e~^/('^+^) (see Fig. 3) is not due to the 
presence of an eigenvalue nearby, for the deformed operator. 

Proof of Proposition [TT9t We shall prove the lower bound (|1.17p in two steps. We first observe 
that, in view of Proposition 14. H it is sufficient to show that there exist positive constants ci, C2 
such that 

a{H,) n {z G C; ciRez < |Imz| < = , (5.1) 

for e > small enough. Indeed, assume for the moment that (j5.ip holds. If A € M is such that 
|A| > C2 e~^, the proof of Proposition 14.11 shows that K{e, A) = \\{H^ — iA)"^|| < Ce^/^ for some 
C > 0, hence any z £ o'(i^e) with jlmzj > C2€~^ necessarily satisfies Rez > C~^€^^^'^. On 
the other hand, the spectrum of H^^ in the strip |Im2:| < C2 is concentrated in the sector 
I Imz| < ci Rez by (jS.ip . hence to prove the lower bound (jl.lZp it is sufficient to consider the 
smaller strip |Imz| < e~'^ . We thus take A £ R with |A| < e^*^ and we estimate K(e, A) as in 
the proof of Proposition 14.11 Using exactly the same notations, we find for all j G N: 

C,(6,A) > ^(2'^~^ + -L- inij2'^f{2^x)-2'^eX\) , (5.2) 

see e.g. (j4.22p . As usual it is sufficient to consider large values of j, in which case 2^^f{2^x) 
\x\~^ for all X e Kj = {x e R; \ < \x\ < 1}. If 2'^-'e|A| < 1/4, the infimum in (15. 2p is larger 
than 1/2 for large j, hence 
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If 2'=Je|A| > 1/4 and |A| < e"'^', then 2^=^ > Ce""'-^ hence Cj(e,A) > 22^-5 > Ce''^'. Thus we 
have shown that K(e, A)^"*^ > C valj Cj{e,X) > Ce"'' if |A| < , from which we deduce that 
any z G cr(H^) with |Imz| < e~'^ necessarily satisfies Rez > C^^e~'^ . Summarizing, we have 
shown that there exists M5 > such that has no spectrum in the half-plane {Re z < M^e"'^ } 
if e is sufficiently small, which is the desired result. 

It remains to prove (15. ip . Clearly this cannot be done using only the pseudospectral estimates 
of Proposition HTH because by Corollarv l4.2l the domain {z G C ; ci Rez < | Imz| < C2 e~^} meets 
the pseudospectrum of as e ^ 0. Instead we shall implement a complex deformation method 
using the group of dilations {Ug(p){x) = e^^'^(f{e^x), which are unitary operators when G M. 
If / is given by (I1.16|] . the multiplication operator by {i/e)f{x) is a dilation analytic perturbation 
of the harmonic oscillator hamiltonian Hoo = —d'^ + x"^, according to Definition 8.1 in [2]. This 
implies that 

defines an analytic family of type (A) in the strip S = {9 e C ; \ lm6\ < tt/4} (see [H]), with 
common domain D = D{Hao)- In particular the spectrum of H^(6), which is always discrete, 
does not depend on 6* G 5". 

We choose 6 = itk, where tk = observe that the operator H^{itk) is still 

maximal accretive. Indeed, for all (x,^) G M^, 



^2 + e^'^x^x^ G |z G C ; -2tk < arg(z) < 2tfc| C {z G C ; Re(z) > 0} , 

^-ktk< arg(z) c {zeC; Re(z) > 0} , 



i (1 + e-2^*fex2)'=/2 

G <! z G 



e |1 + e^^ifex^l'^ 

hence Re(if(itfc)c/? , i^) > for all ip £ D. It is therefore sufficient to reproduce part of the 
analysis done in Proposition 14.11 to estimate ||(-f^e — iA)~^|| for A G M. Consider a partition 
of unity Xo + xL = ^ with suppxo C (-1,1) and xo = 1 on [-1/2,1/2]. Applying the IMS 
localization formula ()4.4p to the operator Q = e^***^ (ff^ (it^ ) — iA), and using the fact that 

Wd^uf < ^Re{{H,{itk)-iX)u,u) < V2\\u\\\\{H,{itk) - iX)u\\ , 

cos(2ifcj 

we obtain 

3\\{H,{itk) - iX)uf + C\\uf > \\{H,{itk) - iX)xouf + m,{itk) - iX)xoouf . (5.3) 

We need to estimate both terms in the right-hand side. 

i) The localization of uq = xou in the interval (—1,1) implies 

||uo|| \\{H^{itk) - iX)uo\\ > Im{{H^{itk) - iX)uo , uq) 

> ^ Re ((/(e'^'^x) - Ae)no , uq) - sm{2tk)\\d,uof 

> - inf (Re f{e'^>=x) - eX)\\uof - tM'^tk)\\uo\\\\iH,{itk) - iX)uo\\ . 

e |a:|<l V / 

Since |1 + e^**'^^^! < 2 when < 1 and —kt^ < arg(/(e***x)) < 0, we have 

Re /(e**^x) > > ^ , for |x| < 1 . 
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If we assume that e|A| < 02'= 2 2 , we thus arrive at the lower bound 

2\\{H,{itk) - iX)xou\\ > (1 +tan{2tk))\\{H,{itk) - iX)xou\\ > ce-^Uxo^^ll , 
where c = — 2~~^ > 0. 

ii) On the other hand, setting u^o = XooU and assuming that Uoo 7^ 0, we observe that 
\\{H^{itk) - iX)uc 



^00 II ^ 



{H^{itk)u 

00 ) ^00 / • > 
— — lA 



I ^00 1 1 ^ 



> inf I2: — zA| , 



ll'^oo II 

where is any sector in the complex plane that is large enough to contain the quantity 



l^^2it,^2\k 

for any Mqo £ ^ with suppuoo C M \ [— ^, If we denote 

Si, = iarg(e2**'= + ^) G (0,tfc) , so that k6k > ^ max^arg((e^'*'= + a;^)''/^) , 
we have f > f — k{tk — S^) > f — kti, > 2tfc, hence we can choose 

Sk = £ C; -2tk < argz < ^ - k{tk - Sk)^ , 

and there exists Cfc > such that inf^gg^ \z — iX\ > Ck\X\. 

If we now combine the lower bounds of both terms in the right-hand side of (15.30 . and if 
we use in addition the fact that \\{H^{itk) — iX)u\\ > cos(2tfc)||u||, we conclude that there exists 
ci > such that 

\\{H^{itk) - iX)u\\ > C minj^ , Cfe|A||||n|| > — ||n|| , for all u G P , (5.4) 

when A G M satisfies |A| < C2 e~^. Thus the deformed operator H^iitk) is sectorial in a neighbor- 
hood of the origin in the complex plane, uniformly for all e E (0, 1]. In particular, if z = -|- iA 
with < ciji < |A| < C2 e~^, it follows from (j5.4p that 

\un(; ^_ ^ \\[HStk)-iXr^\\ , , 

hence z ^ a{H^{itk)) = a[H^). This proves (jS.ip . since we already know that has no spectrum 
in the half-plane Re z < 0. □ 



Although Propositions 14.11 and 11.91 provide the necessary information for comparing the 
spectral bound S(e) and the pseudospectral quantity ^'(e), the analysis can be pushed further 
in the particular example (11.160 . As it is the case for ^(e), the actual behavior of S(e) as e ^ 
results from a competition between contributions due to the (unique) critical point of /, and 
to the behavior of f{x) as \x\ — > 00. Instead of giving a complete proof, we briefly sketch the 
main idea, which is adapted from the usual semiclassical techniques of harmonic approximation 
(see for example [9.]) and complex deformation (see for example |181 121j). We also verify our 
conclusions using numerical simulations. 

We first write eH^ = -h?dl + V{x,h?) where h = e^^"^ and V{x,h^) = h^x"^ + if{x), and 
we remark that eH^ is not exactly a semiclassical operator, because the potential V depends 
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on /i in a nontrivial way which makes the comparison with quadratic approximations at critical 
points more difficult. In the sector {z £ C; |argz| < vr/4}, the function V{z,e) has exactly 
three critical points: z = 0, and z = iz^, where 



'ik\'i^ -11/2 /ik\^ , 1 



The quadratic approximations of are respectively: 



i) At z = 0: F° = -dl + ^ + (1 + '-^^)x^- The eigenvalues of are 

/^n(e) = - + (2n + , n e N , where = (l + 7^/"(0)^ 

ii) At z = ±z,: Hf = -dl + A + ^^?2;2, with 

a z , , A: + 2 , 2 k + 2/ik\'ii' 

p. = , ,n -/(..) = —..= + ^ = —(-) -1, 



0? = l + 5jr(--.) = ('= + 2)^ = (t + 2)(l-(^ 

The eigenvalues of Hf^ are f^(e) = + (2n + n G N. 

After a complex rotation, these quadratic approximations are transformed into real harmonic 
oscillators and H^, respectively. The actual complex deformation has to be chosen carefully 
so that the non-self-adjoint corrections to (resp. H^) still allow for a good control of the 
resolvent norm of the deformed hamiltonian when z is close to ^5^(e) (resp. f^(e)). This 
requires an accurate analysis, which is left for a future work. In any case, on the basis of the 
arguments above, we expect that the full operator has two sets of eigenvalues {/Lt„(e)} and 
{t'„(e)} which satisfy: 

M„(e) = ^ + (2n + + 0{l) , n = 0, 1, 2, . . . (5.5) 
vn{e) = D, + {2n + l)a, + 0{e'') , n = 0,1,2,... (5.6) 

This leads to the following conjecture: 

Conjecture 5.1 Fix k > and let f : M ^ M be as in ()1.16p . Then the spectral bound of 
satisfies, when e > is small, 

E(e) = min{Re(Mo(e)) , MM^))} , (5-7) 
where Ho{e) is given by (j5.5p and z^o(e) by (|5.6p . 

Conjecture 15.11 predicts the following asymptotic expansions as e ^ 0: 

r i//c\i/2 



2 Ve 



+ 0(1) , if A: < 2 



k + 2 / k\-2'^ 



) cos(7rzy) + V/c + 2 - 1 + Oie") if k>2 



k \2eJ 

with V = [k + 2)~^, in full agreement with Proposition 11.91 
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To test these predictions, we have computed an approximation of S(e) using a standard 
finite difference scheme, for three different choices of k and for values of e ranging from to 
2~^*. Calculations were made on a large interval [—L,L] with Dirichlet boundary conditions, 
using at least 500 grid points. The results are collected in Fig. 4, where the dots represent the 
values computed numerically and the solid lines the theoretical predictions, namely Re(/XQ(e)) 
for k < 2 and Re(z^Q(e)) for k > 2. The agreement is excellent and leaves little doubt on the 
validity of Conjecture 15. 1[ 



6 







k = 1 




k = 2 ^"^"^-^ 




= 4 






o ^'^"--^ 
o 



o 



-12 -10 -8 ^ , , -6 -4 -2 

log(e) 



Fig. 4: The logarithm of S(e) is represented as a function of log(e) for three choices of the parameter k. The 
dots correspond to the values computed numerically, and the solid lines to the theoretical predictions given by 
Re(At!3(e)) {k = 1,2) and Re(i/g(e)) [k = 4). To improve readability, the upper curve has been shifted up by two 
units and the lower curve shifted down by two units. 

The heuristic arguments above give a correct prediction not only for S(e), but also for all 
eigenvalues of with sufficiently low real part. Let us first consider the case when /c = 4 and 
e = 2~^^. In Table 1 below, we have listed the five eigenvalues of with smallest real part, 
together with the approximate eigenvalues ^^^(e) and f^(e) for n = 0, . . . ,4. It is obvious that 
f^(e) is an excellent approximation of A„(e). 



n 


A„(e) 











106.18 60.48 i 


512 2.616- 


10^ i 


106.18 60.48 i 


1 


111.02 60.51 i 


1536 + 2.606 


• lO^i 


111.06 60.50 i 


2 


115.83 60.54 i 


2560 + 2.596 


• 10^ i 


115.93 60.5H 


3 


120.61 60.58 i 


3584 + 2.586 


• 10^2 


120.80 60.53 i 


4 


125.36 60.63 i 


4608 + 2.575 


• lO^i 


125.67 60.54 i 



Table 1: The first few eigenvalues of He for k — A and e = 2 together with the predictions and Vn{e). 

We next consider the case when k = 2, which is especially interesting because both Re{fi^{e)) 
and Re(i/°(e)) grow with the same rate 0(e~^/^) as e ^ 0. However we see from (|5.5p . (|5.6p 
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that Re(/x°(e)) ^ (2n + l)/(2e)i/2 while Re{u^{e)) « (2/e)V2 + 4n + 1. It follows that the 
eigenvalue of with minimal real part is well approximated by ^o(e), and thus corresponds to 
a semiclassical mode near the origin, while the subsequent eigenvalues are well approximated 
by fni^) — 0, 1,2, . . . and therefore correspond to semiclassical modes near infinity. This 

prediction is clearly confirmed by the values listed in Table 2 below. It is also possible to 
verify that the eigenfunctions corresponding to A„(e) are nicely approximated either by an 
eigenfunction of the quadratic approximation at z = (if n = 0) , or by an eigenfunction of 
the quadratic approximation Hf^ at z = ztz^ (if n = 1, 2, . . . ). 



n 


An(e) 









44.54 + 4051 i 


45.26 + 4050 i 


91.50 + 90.52 i 


1 


91.50 + 90.52 i 


135.78 + 3960? 


95.48 + 90.54 z 


2 


95.47 + 90.54 i 


226.30 + 3869 i 


99.45 + 90.57 i 


3 


99.48 + 90.56 i 


316.82 + 3779 i 


103.43 + 90.59 i 


4 


103.2 + 90.58 i 


407.34 + 3688 i 


107.41 + 90.6H 



Table 2: The first few eigenvalues of for k — 2 and e = 2 together with the predictions fiZ{^) and !^JJ(e). 

Finally we turn our attention to the case when k = 1. Here we expect that the eigenvalues 
of with lowest real part are well approximated by /x^(e), for n = 0, 1, 2, . . . . This is indeed 
the case, although the spectrum of seems more difficult to compute in this regime, due to 
numerical instabilities. Nevertheless, for e = 2~^^, we were able to check that at least the first 
three eigenvalues of are correctly approximated by //°(e), see Table 3 below. The errors are 
larger than in the preceding cases, but we should take into account the fact that the remainder 
term in ([53]) is 0(1). 



n 


An(e) 









31.46 + 4064 i 


32.00 + 4064 i 


242.63 + 419.00 i 


1 


93.34 + 4002 i 


96.02 + 4000 i 


246.09 + 419.01 i 


2 


153.2 + 3940 i 


160.0 + 3936 i 


249.55 + 419.01 i 



Table 3: The first few eigenvalues of He for fe = 1 and e = 2 together with the predictions Mn(e) and 

A Appendix 

This appendix gathers the proofs of Lemmas 11.11 11.31 11.71 and provides a variant of Lemma ll.7l 
which is needed in Section 13.31 

We first give a proof Lemma II. 1] which relates the quantities S and ^ to the norm of the 
semigroup e^*"^ for a sectorial operator, as an application of the Laplace transformation. 

Proof of Lemma 11.11 Let ^ be a maximal accretive operator in a Hilbert space X, with 
numerical range Q{A) contained in the sector = {z £ C; |arg2:| < ^ — 2a} for some 
aG(0,f]. 

i) Assume first that He"*"^!! < Ce~^* for all t >0. Then the resolvent 

/•oo 

(A-z)-'^ = / e-*^e*Mt 
Jo 
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is defined (at least) for all z G C with Re z < fi, hence S > /x. Moreover, taking z = iX ^ il 
and using the fact that e"*"^ is a semigroup of contractions in X, we find 



WiA-tX) 



< 



foo poo 

/ ||e-*^||dt < / min{l, Ce-'^^jdi 
Jo Jo 



1 + log(C) 



Taking the supremum over A € M, we conclude that ^' > /u/(l + log(C)). 

ii) Conversely, to estimate the semigroup e"*"^ in terms of S or ^, we use the inverse Laplace 
formula 

e-*^ = — [ {A - z)-^ e"' dz , 
2vrz J^i^l,a) 

where < /i < S and r(^, a) = r_ {fi, a) U ro(^, a) U r+(^, a) is the polygonal contour defined 
by 

To{fi, a) = 1^ S C ; Re z = /i , I arg 2:| < — — a| , 
T±{^,a) = ja: S C ; Re 2 > /i , arg z = lb ^— — I , 

and oriented from Imz = — oo to Imz = +oo. Note that r(/i, a) lies entirely in the resolvent 
set of A by construction. Since Rez = fi when z G ro(/U, a), we easily obtain 



/ {A-z)-'^e-''dz <N{A,fi)^^ 
Jro{fi,a) tana 



-fit 



where N{A,fi) = sup{||(A — 2;)^"'^|| ; Re(z) = /u}. On the other hand, any z G r+(/i,a) can be 
parametrized as z = x + {ix/ tana) with x > fi, and the resolvent at this point can be estimated 
as follows: 



\\iA-z)-'\\ < 



1 



< 



We thus find 



(A-z) 



-1 -zt 



dist(z,e(A)) ~ dist(z,Ao) x ' 

1 e-t"^ 



dz 



< 



-tx 



dx 



< 



tan a tan a fit 



and the contribution of r_(/i, a) is estimated in exactly the same way. Collecting these bounds 
are using the fact that He"*"^!! < 1 < 1/ tana, we arrive at 



-tA\ 



< 



1 



Trtana 



nN(A,n)e-^'^ +mm(TT,^}] < — - — ( fiN (A, n) + 27r) e' 
L fJ.t ) J vr tan a V / 



ill) Assume that < < ^. Using the second resolvent formula and the definition of we 
find for all A G M: 



\\{A-fi-iX)-'\\ < 



< 



1 



ll(^-^A)-^ll 
I - fi\\{A - iX)-^\ ^ ^--/z 



hence N{A,fi) < - 



□ 



We next give a proof of Lemma 11.31 which illustrates the pseudospectral nature of the 
quantity ^'(e). 



Proof of Lemma 11.31 
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i) If Re(z) < 0, we know that \\{H^ - z)~^\\ < 1/ dist(z, e(F,)) < 1. If < Re{z) < K^{e) for 
some K G (0, 1), we have by Lemma 1 1.1 1 - iii): 



^'(e) - Re{z) - ^'(e) 1-k " 1-k 

Thus \\{H^ — is uniformly bounded by (1 — /t)~^ when Ke{z) < K^'(e). 

ii) We argue by contraposition. Fix Kq > 1, G N, and assume that \\{H^ — z)~^\\ < KQe~^ 

whenever Re{z) < fi^. We shall show that there exists K > 1 (independent of e) such that 

< K^{e) (l + log ^-(6) + log(e-i)) , for ah e G (0, 1] . (A.l) 

Note that (lA.ip is automatically satisfied if /i^ < 1, hence we assume from now on that //e > 1. 
Applying Lemma [TTTI -ii) with fi = fj.^, and using the fact that a = 0{e) when A = H^, we see 
that ||e-*^^|| < C{H^,|l^)e-^'^* for ah t > 0, where 



Trtana 

for some Ki > 1. Next, using the lower bound on ^ given by Lemma ll.H -i). we obtain 

/i. < ^'(e)(l + logC(F„/x,)) < vI/(e)(^l + log(^,)+log(A'ie-^-i) 

The desired bound (jA.ip is now a direct consequence of the following elementary result: 
Claim: // ^ > 1, ^' > 1, and C >1 satisfy fj, < ^{C + log fi), then 

/i < ^{C + 2 log 2 + 2 log ^' + log C) . (A.2) 

Indeed, the hypothesis implies that /U < /Uq, where ^ 1 is uniquely determined by the relation 
Ho = ^{C + log/Uo). Since log no < y/JIo-, we have //q ^ ^(C + ^/JIo), hence 

VW) <\{^ + \/^2^^4C¥) < ^ + Vc^ . 

Using this bound and the fact that /i < '^{C + log/xo), we easily obtain ()A.2p . This concludes 
the proof of the claim, hence of Lemma 11.31 □ 

We now briefly recall the comparison argument with semiclassical models which yields the 
lower bound (jl.lip . 

Proof of Lemma 11.71 Let V{x^e) = + ^ ^ . By the Min-Max principle, it is sufficent to 
find a constant M3 > 1 such that 

{H,u, u) = I (\dM^ + V{x,e)\u\^\ dx > -Ml^^ ^ for ah u G P , (A.3) 
7m ^ ^ e^'' 



and for any e G (0, 1] a nonzero function (p^ G C^(M) C V such that 

{H.^e,^e) < ^\Wef • (A.4) 

Hypothesis 11.61 ensures that there exists L > such that 



f{xf > , for |x| > L 

2 xpC'+i) ~ 
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and that / has a finite set {xi, . . . ,xn} of (nondegenerate) critical points. Hence there exists a 
partition of unity Ylf=o = ^ such that Xj ^ C'^(M) for j G {1, . . . , N}, suppxo C (—00, — L)U 
(L,+oo), and such that / has exactly one critical point (namely, Xj) in the support of Xj for 
j G {1, . . . , N}. The IMS localization formula of [2] provides a constant C > such that 

/ + y(x,e)|u|2) dx > y / (id^Ujl"^ + V{x,e)\uj\'^) dx - C\\uf , (A.5) 

with Uj = XjU- By construction, for j £ {1, . . . , N}, there exists Cj > such that 

V^(2;,e) > — ^ — > K , for all x G suppxi • 

Hence the comparison with the harmonic oscillator hamiltonian —9^ + c|e^^(x — xj)'^ implies 
that {HeUj , Uj) > Cj e~^\\uj\\'^ for all j G {1, . . . , A^}. For j = 0, we notice that 

^2 2_ 

V{x, e) > x^+ ^^2\x^2(^k+i) - ^0 ^ ' ^ ^ suppxo C M \ L] . 

Thus the first term {H^uq,uq) is bounded from below by c^e^'^. Summing up all the lower 
bounds of the terms in (jA.Sp . and recalling that < {H^Uj ,Uj), we obtain ()A.3P provided 
M3"^(l + C) < min{cj ; < j < A^}. For the upper bound (jA.4p . it suffices to take (feix) = 
ip{e''x) with if G C^{{1, 2)) and ip^O. □ 

Following exactly the same lines, one can also prove the following result, which is needed in 
Section I3.3t 



Lemma A.l Let (3 : M ^ M-|_ be the function defined in Section \3.3\ and depicted in Fig. 1. 
There is a constant Mq > such that, for any u £ T> and any e G (0, 1], 

j {\do,u\^ + x^\u\^ + ^ f {xf\u\^^ dx > ^\\u\\l2 , where D = -^-^ • (A.6) 

Proof: We assume that the parameter ^ > entering the definition (j3.18p of the function [3 is 
large enough so that all critical points of / are contained in the interval [—A + 1, A — 1], and so 
that 

f'{xf > 2^rc^2(k+i) ' whenever \x\ > A . (A.7) 

We now consider the potential V{x, e) = x^ + /^(^)/ (^) g^j^j introduce like in the proof 
of Lemma [1.71 a partition of unity Ylf=o^'j = ^ with Xj ^ C'o°(R) for all j G {!,..., A^}, 
suppxo C {—00, —A) U (^,+00), and such that / has exactly one critical point in suppxj for 
j G {!,..., A^}. Again the IMS localization formula provides a constant C > such that 




with Uj = Xj'ti- For j G {1, . . . , A^} we find as in the proof of Lemma 11.71 a constant Cj > such 
that 

J^(\d,,Uj\^ + V{x,e)\ujf)dx > ^\\ujf . (A.8) 
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For j = 0, we use the definition (j3.18p and the condition ()A.7p to find an appropriate lower 
bound on V{x, e). If cq > is sufficiently small, then for any x £ suppxo we have either 

A<\x\<B,, so that V(x, e) > x^ + ,° ^ > cq e"^/^ , 

or 

k 

n Bne~'^k'^ ?_ 

\x\ > , so that T^(x,e) > x'' + , ,,,, , > cq e ^+4 . 

Since = < ^, it follows that V{x, e) > cq if x G suppxo, hence (lA.Sh holds for j = 
too if we replace the exponent 1/2 by u. We conclude as before. □ 



References 

[1] J. Aguilar, J.M. Combes. A class of analytic perturbations for one-body Schrodinger Hamil- 
tonians. Comm. Math. Phys. 22 (1971), 269-279. 

[2] H.L Cycon, R.G Froese, W. Kirsch, and B. Simon. Schrodinger operators with application 
to quantum mechanics and global geometry. Text and Monographs in Physics, Springer, 
1987. 

[3] E.B. Davies. Non-self-adjoint differential operators. Bull. London Math. Soc. 34 (2002), 
513-532. 

[4] N. Dencker, J. Sjostrand, and M. Zworski. Pseudospectra of semi-classical (pseudo)- 
differential operators. Comm. Pure Appl. Math. 57 (2004), 384-415. 

[5] L. Desvillettes and C. Villani. On the trend to global equilibrium for spatially inhomoge- 
neous kinetic systems: the Boltzmann equation. Invent. Math. 159 (2005), 245-316. 

[6] Th. Gallay and C.E. Wayne. Invariant manifolds and the long-time asymptotics of the 
Navier-Stokes and vorticity equations on M?. Arch. Ration. Mech. Anal. 163 (2002), 209- 
258. 

[7] Th. Gallay and C.E. Wayne. Global stability of vortex solutions of the two-dimensional 
Navier-Stokes equation. Comm. Math. Phys. 255 (2005), 97-129. 

[8] M. Hager and J. Sjostrand. Eigenvalue asymptotics for randomly perturbed non-selfadjoint 
operators. Preprint. 

[9] B. Helffer. Semi- classical analysis for the Schrodinger operator and applications. Lecture 
Notes in Mathematics 1336, Springer, Berlin, 1988. 

[10] B. Helffer and F. Nier. Hypoelliptic estimates and spectral theory for Fokker-Planck Opera- 
tors and Witten Laplacians. Lect. Notes in Math. 1862, Springer, 2005. 

[11] B. Helffer and J. Nourrigat. Hypoellipticite maximale pour des operateurs polynomes de 
champs de vecteur. Progress in Mathematics 58, Birkhauser, 1985. 

[12] B. Helffer and J. Sjostrand. Resonances en limite semiclassique. Mem. Soc. Math. France 
(N.S.) No. 24-25 (1986). 

[13] F. Herau and F. Nier. Isotropic hypoellipticity and trend to the equilibrium for the Fokker- 
Planck equation with high degree potential. Arch. Ration. Mech. Anal. 171 (2004), 151-218. 

[14] F. Herau, C. Stolk and J. Sjostrand. Semiclassical analysis for the Kramers-Fokker-Planck 
equation. Comm. Partial Differential Equations 30 (2005), 689-760. 

[15] P.D. Hislop and I.M. Sigal. Introduction to spectral theory. With applications to Schrodinger 
operators. Applied Mathematical Sciences 113, Springer, 1996. 



37 



[16] L. Hormander. The analysis of linear partial differential operators. Springer, 1985. 

[17] T. Kato. Perturbation theory for linear operators. Grundlehren der mathematischen Wis- 
senschaften 132, Springer, 1966. 

[18] A. Martinez. An introduction to semiclassical and microlocal analysis. Universitext, 
Springer, 2002. 

[19] K. Pravda-Starov. A general result about pseudospectrum for Schrodinger operators. Proc. 
R. Soc. Lond. A 460 (2004), 471-477. 

[20] L.P. Rothschild and E.M. Stein. Hypoelliptic differential operators and nilpotent groups. 
Acta Mathematica 137 (1977), 248-315. 

[21] J. Sjostrand. Singularites analytiques microlocales. Asterisque 95, 1-166, Soc. Math. France, 
1982. 

[22] J. Sjostrand. Some results on non-self-adjoint operators, a survey. Preprint available on 
|http://arxiv.org/pdf/0804.3701vlj 

[23] L. Trefethen. Pseudospectra of linear operators. SIAM Rev. 39 (1997), 383-406. 

[24] C. Villani, Hypocoercivity, Memoirs of the AMS, to appear. 

[25] C. Villani. Hypocoercive diffusion operators. International Congress of Mathematicians. 
Vol. Ill, 473-498, Eur. Math. Soc, Zurich, 2006. 



38 



